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Abstract 


A new software package has been developed for use with Micro and 
Mlnlcoaiputers, particularly Digital Equipment Corporation's LSI-ll/PDP-ll 
series. The package has a nueber of Fortran-callable subroutines which 
perform a variety of frequently needed tasks for biomedical applications. All 
routines are well documented, flexible, easy to use and modify, and require 
minimal programmer knowledge of peripheral hardware. The package Is also 
economical of memory and CPU time. A single subroutine call can perform any 
one of the following functions: (1) Plot an array of Integer values from 

sampled A/D data; (2) Plot an array of Y values versus an array of X values; 
(3) Draw horizontal and/or vertical grid lines of selectable type; (4) 
Annotate grid lines with user units; (5) Get coordinates of user-controlled 
crosshairs from the terminal for Interactive graphics; (6) Sample any analog 
channel with prcqr.» selectable gain; (7) Walt a specified time Interval; and 
(8) Perform random access I/O of one or more blocks of a sequential disk file. 
Several miscellaneous functions are also provided. These routines are modular 
and easily changed, and are especially applicable for uses In biomedical 
research laboratories such as NASA's where adaptability Is Important and soft- 
ware development time Is limited. Complete source code listings, example main 
programs, and sample output are Included. 
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Introduction 


Numerous experiments In biomedical research laboratories Involve the 
collection and display of sampled analog data using a mini or microcomputer. 
Often, however, many experiments which Inltlaly seem simple, turn out to 
Involve a considerable Investment In programing effort because of the 
difficulty In finding suitable graphics packages, or In using such devices as 
analog to digital (A/0) converters In a high level language, such as Fortran. 
Most data acquisition routines (except for those Involving realtlvely low 
sampling rates) tend to be application specific, and are written In assembly 
language by Individuals who must have a fairly Intimate knowledge both of the 
particular A/D systems used and of the assembly language of the computer. In 
addition, biomedical data often seems particularly well suited for graphical 
displays, because of the complexity of physiological data and also because of 
the ease of Interpretation of graphs by medical or other personnel without an 
extensive computer or statistical background. Graphics software packages are 
available from various sources, but most are not particularly well suited for 
use with microcomputers. Many of the commercially available packages are 
quite large and have features that are much more extensive than needed In many 

applications. 

Some of the features are used rarely. If at all, but their presence still 
requires a rather large amount of main memory and disk space In the system. 
The architecture of some of these packages can be quite complicated, and with 
many of their routines rather poorly documented. It can be a formidable task 
for a programmer to remove the extra, unneeded features from such a package 
without interfering with the basic functions that are required. As a result, 
these large packages are often only suited for larger computer facilities, 
where large amounts of memory are available and where numerous users may make 
use of the myriad of features offered. 

Another problem with some graphics packages Is that they are not 
particularly easy to use with real time Interactive programs. A1 though fine 
for offline plotting appllatlons, they may be too slow or too difficult to use 
In cases where real time data acquisition and control of an experiment are 
taking place. In addition, some packages are not sufficiently well 
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doc unen ted, particularly with comments In the subroutines or with examples, to 
allow them to be used or modified easily by programmers without extensive 
knowledge of the hardware being used. Sometimes, It can take a naive 
prog runner several weeks to write one fairly simple program, because of the 
time required to discover the pecularltles of analog to digital converters and 
graphics terminals or plotters. The plethora of functions available In some 
packages also seems only to cause confusion among some users, since many of 
the functions are redundant and do not have to be used at all. 

Because of these problems with commercially available software, a new 
data acquisition graphics package was developed for use at the NASA Johnson 
Space Center Life Sciences Laboratories. It was designed for use specifically 
with micro or minicomputers, particularly Digital Equipment Corporation' s 
LSI-11 and PDP-11 series, using the RT-11 operating system. Desired features 
of the package Include the following: (1) It should be relatively small In 

size so that It can be used easily with microcomputers having limited memory 
(56k bytes or less). (2) It should Include the most frequently needed 
graphics, analog to digital conversion (A/D), and miscellaneous capabilities. 
However, extra functions which are not commonly needed or which can be 
performed with alternate methods, and which require too much memory , should 
not be included. (3) The package should be modular, flexible, easy to learn, 

use, and modify. Additionally, use of the package should require only a 
minimal knowledge of the peripheral hardware. The package should be adaptable 

to various types of Tektronix terminals or possibly to X-Y plotters. It 
should also be useable with Fortran main programs. (4) In order to meet the 

previous objective, the package should be written primarily In a high level 
language, and should be well commented and otherwise documented as well. (5) 

Finally, It also should be able to handle moderately fast analog to digital 
conversion and displays for real-time applications. 

With the above objectives In mind, a software package was developed which 
has the following characteristics. First, It consists of a nunber of Fortran- 
callable subroutines which perform all necessary tasks required for 
Interfacing with an A/D system or various types of Tektronix terminals. Only 
a few of these routines would need to be modified If a different type of 
terlnal or X-Y plotter was used. Further, most of the subroutines themselves 
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arc written In Fortran and all are fully consented so that they can be easily 
Modified by others. If necessary. In addition, a few other staple 

Miscellaneous routines were written to perfona other tasks used frequently In 
laboratory applications. A list of all of the routines In the package Is 
shown In Table 1. 


Routine 
GRINIT - 
MPLOT 

COPY 

ERASE 

CHRSIZ - 
GINPUT - 

ARYPLT - 

XYPLOT - 


GRIO 

AMOTAT - 
BELL 

WAIT 

ISAMPA - 
DISKIO - 


TABLE 1 

SUBROUTINES IN THE SOFTWARE PACKAGE 
Description 

Initialize graphics paraMeters for other routines 
Change position (move cursor), go to alphanuwrlc Mode, or 
draw vector 

Make hard copy of Tektronix terminal screen A erase It If 
desl red 

Erase Tektronix screen without copying 
Change character size 

Switch to graphic Input Mode, display A get coordinates of 
user-controlled crosshairs 

Plot an Integer array of Y values (A/D saMples, etc), with 
straight lines connecting points, with variable scaling 
Plot an array of real Y values vs an array of real X 
values, with straight Interconnecting lines, with variable 
scaling 

Draw grid lines (selectable type) over desired area 
Label plot axes with user units at some or all grid lines 
Ring terminal bell /beep (Variable duration A modulation 
control ) 

Walt a desired period of ttae (1/60 sec resolution) 

Sample an A/D channel with selectable gain A channel master 
Perform random access multiblock binary I/O to a sequential 
disk file 
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A few of these routines tlso use subroutines from the Fortran library In 
the system to access the system line frequency clock, to send characters to 
the terminal through the terminal handler, to perform double word Integer 
arithmetic, and do disk I/O. These few routines In the package would have to 
be modified If a different computer or operating system was used. In 

addition, routines MPLOT , COPY, ERASE, and possibly GRINIT would have to be 
changed If a different terminal or plotter was used. However, such changes 
should be fairly straightforward once the characteristics of any particular 
terminal were determined. 

The package will presently support the following hardware on an LSI-11 
series microcomputer. First, analog to digital conversion can be performed, 
with LSI-11 compatible A/D systems manufactured by ADAC, Data Translation, or 
DEC, using different modes. Programmable gain and random channel selection 
are supported. A throughput rate of 4000 samples per second or more can be 
accomplished If little or no computation Is performed between samples. 
Furthermore, all of the Tektronix 4000 series graphics terminals can be used 
with the package. Features of some of these terminals which are supported, 
but not required, Include: (1) 12-bit addressing for X and Y coordinates, (2) 

variable character sizes and dot-dash line types (using the hardware built 
Into the terminals), (3) graphic tnp«»t mode for Interactive applications (so 
that a program can determine the positions of user controlled crosshairs), and 
(4) write-through mode for non-stored (refreshed) displays. Additionally, 
several different types of X-Y plotters could be used with the package, 
although routine MPLOT would have to be modified to accommodate the 
requirements of the particular plotter used. 

The remainder of this special report consists of the following. First, 
some simple example main programs are given In the next section In order to 
demonstrate the use of, and capabilities of, some of the routines In the 
package. Graphical outputs from the sample programs are also provided. 
Following that section, each subroutine In the package Is discussed 
separately, with Information on calling conventions, parameters (arguments) to 
be passed, and any restrictions or additional details concerning Its 
operation. Source listings of each subroutine are also given In each section. 
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Additional examples and other programs are provided In the references 
below. 
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INSTALLATION 


All of the subroutines In this package are callable by Fortran program 
and are designed to be used with DEC'S RT-11 operating systeai. For conven- 
ience, some or all of then My be placed In a library on the sane or a 
different disk as the System Subroutine Library (SYSLIB) . 

For example, In order to create a library called LABLIB containing 
routines GRINIT, MPLQT, ARYPLT, and ISAMPA, do the following. First, put the 
source code for each routine In separate files on the default device DK:, each 
with the name of the subroutine: GRINIT. FOR, MPLOT.FOR, ARYPLT. FOR, and 

ISAMPA. MAC. Second, compile and assemble these routines with the RT-11 
(Version 4) Fortran compiler (Version 02.1 or later) and Macro assembler 
(Version 4 or later): 

FORT/WARN/LIST:DK: (GRINIT, tf>L0T, ARYPLT) 

MACRO /CR/L 1ST: DK: ISAMPA 

Next, create the actual library file LABLIB. OBJ from the Individual 
object files: 

LIBR/CREATE LABLIB 

Files? GRINIT.MPLOT, ARYPLT, ISAMPA 

The same procedure may also be followed for more subroutines, but you 
should Include only six or fewer file names on any command line. If you want 
to put routines from more than six files In the library, then Include the 
"/PROMPT" option after the command "LIBR", then specify six or fewer file 
names per line. After the last line, type two slashes (//) to terminate the 
file name entry. Refer to the RT-11 System User's Guide for more Information. 

You will note. In the descriptions of several of the subroutines, that 
they In turn call other subroutines. For Instance, MPLOT and WAIT are called 
by a number of other routines. Make sure that all of the subroutines which 
are needed (either directly or Indirectly) are Included In the library which 
you build. Otherwise, an undefined global error will occur when you attempt 
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to link the programs. In addition, several DEC suppllad routines from the 
System Subroutine Library are used, and these must be present In SYSLIB. Such 
routines Include ITTOUR (called by MPLOT, COPY. ERASE, CHRSIZ, BELL, etc.) and 
GTIM and the Integer*4 routines (called by WAIT). Refer to the descriptions 
of the appropriate routines In this report for more Information. 


Sample Main Programs 


This section contains several simple programs which demonstrate some of 
the capabilities of the subroutines In this package. By referring to the 
program listings and sample output which follow, and running these programs 
yourself, you should be able to learn how to use the subroutines to perform 
other similar functions. 

The first example program Is AOTEST. This uses only the routine ISAMPA 
from the subroutine package. In addition, ITTINR and IPOKE from the System 
Subroutine Library are also used. More Information on these last two routines 
Is In the DEC RT-ll Version 4 Programmer's Reference Manual. 

Program AOTEST does tht following. First, It asks the user for the types 
of ana1og-to*d1g1ta1 (A/0) and dlgltal-to-analog (D/A) converters present In 
the system. Second, the program asks which analog Input channel should be 
sampled, and the programmable gain code to use. After this, the program 
samples the specified channel with subroutine ISAMPA and sends the converted 
signal back out to the D/A converter(s) In the system where an oscilloscope 
can be used to monitor the D/A output. The sampling Is repeated until the 
user strikes the Return key. The D/A output resembles a sampled (choppy) 
version of the Input analog signal, with a scale factor dependent on the 
progranmable gain, A/D range jumpers, and D/A range jumpers selected. 

Program ADTEST may also be used for diagnostic purposes, since It enables 
one to quickly check both the A/D and D/A converters plus certain CPU and 
memory functions, without being particularly difficult to use. In addition. 
It may be easily modified for other purposes. The program listing and sample 
output follows. 
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FORTRAN 10 


Om&UPL *f 
or pocn ouauv/ 
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0001 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0002 

c 

c 

0003 

c 

0004 

0005 

0006 

0007 

0008 

c • 

0010 

0011 

0012 

0013 

C 

0015 

0016 
0017 

c 

0019 

0020 
0021 
0022 

0023 

0024 


PROGRAM AOTEST 

PURPOSE! TEST A/D S D/A BOARD 

METHOD! SAMPLE SELECTED A/D CHANNEL l TRANSFER SAMPLED 
VALUE TO BOTH D/A'S» AND REPEAT UNTIL 
INTERRUPTED BY USER. USER MAY INPUT 10 HZ 
TRIANGLE WAVE FROn SIGNAL GENERATOR* OR ANY 
OTHER SIGNAL* AND D/A'S SHOULD FOLLOW THE INPUT. 
SCALING MAY BE DIFFERENT* HOWEVER* DEPENDING ON 
PROGRAMMBLE GAIN AND VOLTAGE RANGES SELECTED. 

WRITTEN BY! WILLIAM G. CROSIER 
REVISED! 23 DEC. 1981 

SUBROUTINES REQUIRED! ISAMPA 


INTEGER DAC1 

SET UP DEVICE ADDRESSES 

DAC1 * *176760 ! D/A #1 DATA REGISTER 

TYPE **'A/D S D/A CONVERTER TEST PROGRAM — VERSION 3' 

TYPE *, 'ENTER A/D TYPE! 0 FOR ADAC* OR' 

TYPE **' 1 FOR DATA TRANSLATION' 

ACCEPT # * I ADTYP 

IF ( IADTYP.LT .0 .OR. I ADTYP, GT . 1 ) STOP 'ILLEGAL A/D TYPE' 

TYPE ** 'ENTER D/A TYPE! 0 FOR ADAC. OR' 

TYPE **' 1 FOR DATA TRANSLATION' 

ACCEPT # » IDATYP 

IF (IDATYP. LT.O .OR. IDATYP , GT . 1 ) STOP 'ILLEGAL D/A TYPE' 

TYPE *, 'ENTER NUMBER OF D/A CONVERTERS ON BOARD (1-4)' 
ACCEPT A * NDACS 

IF (NDACS.LT. 1 .OR. NDACS. GT. 4) STOP 'ILLEGAL NO. OF DACS' 

TYPE ** 'VERIFY THAT A/D IS JUMPERED FOR BIPOLAR*' 

TYPE **'2"S COMPLEMENT OPERATION. WITH CSR ADDRESS' 

TYPE *.'0F 176770 OCTAL' 

TYPE Y * ' ALSO VERIFY THAT D/A IS JUMPERED FOR BIPOLAR' 

TYPE Y* 'OPERATION* WITH ABDR OF FIRST D/A OF 176760 OCTAL.' 
TYPE *. 'ENTER -1 FOR INPUT CHANNEL NO. TO STOP.' 


C REQUEST A/D CONTROL PARAMETERS 

C 

0025 5 TYPE 1005 

0026 1005 FORMAT ('OA/D INPUT CHANNEL NO. ?'*$) 

0027 ACCEPT Y* ICHAN 

0028 IF ( I CHAN .LT.O .OR. ICHAN.GT. 31) STOP 

0030 20 TYPE Y* ' PROGRAMMABLE GAIN CODE!' 

0031 IF ( I ADTYP . EQ . 0) TYPE 1020 

0033 1020 FORMAT (' 0 = GAIN OF 10' / ' 1 - GAIN OF 5' / 
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OF POOR QUALITY 
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0034 

0036 

0037 

0038 

0039 

0041 

0042 


0043 

0045 

0046 


0048 

0049 

0050 

0051 

0052 


0053 


1030 


C 


C 

c 

c 

c 

40 

C 

C 

c 


50 

c 
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' 2 = GAIN OF 2' / ' 3 * GAIN OF T) 

IF < IADTYP .EQ. 1‘ TYRE 1030 

FORMAT (' 0 ■ GAIN OF 1 ' / ' 1 * GAIN OF 2' / 

2 * GAIN OF 4' / ' 3 * GAIN OF 8'> 

TYPE * > ' GAIN CODE DESIRED (0-3)?' 


ACCEPT * • IPGNCD 


IF < IPGNCD. LT.O .OR. IPGNCD. GT. 3) GO TO 20 


TYPE * » 'HIT RETURN KEY TO TERMINATE SAMPLING I CHANGE A/D' 
TYPE *,' CHANNEL NO. OR PROGRAMMABLE GAIN.' 


DO A/D CONVERSION & OUTPUT DATA TO BOTH D/A'S 
UNTIL USER STRIKES RETURN KEY ON TERMINAL. 


IF UTTINRO .G£. 0) 60 TO 5 ! HAS RETURN KEY BEEN HIT 

I DATA * ISAMPA ( I CHAN, IPGNCD, IADTYP ) ! SAMPLE A/D CHAN. 

IF DATA TRANSLATION D/A, CONVERT CODING TO OFFSET BINARY 
IF i I DAT YP .EQ. 1) IDATA=IDATA+2043 

IDAC = DAC1 ! D/A CONV . ADDR. 

DO 50 K* 1 , NDACS IFOR EACH D/A, 

CALL IPOKE( IDAC, IDATA) ! OUTPUT SAMPLE TO D/A 

IDAC = IDAC + 2 ! ADDR. OF NEXT D/A 

60 TO 40 

END 
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T ype 

Off set 
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Tyre 

Offset 

Name 

Type 

Offset 

DAC 1 

1*2 

000000 

IADTYP 

1*2 

000002 

ICHAN 

1*2 

000010 

IDAC 

1*2 

OOOOlo 

IDATA 

1*2 

000014 

IDATYP 

1*2 
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1*2 
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1*2 
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The second sample program Is called DEMOGR. This program demonstrates 
the use of GRIHIT, CHRSIZ, ERASE, and HPLOT to plot Interesting patterns on 
the terminal screen. Subroutine ERASE also calls WAIT, and the subroutines 
also use several routines from the System Subroutine Library (such as ITTOUR, 
GTIM, and the Integer*4 functions). 

Program operation proceeds as follows. First, graphics parameters are 
Initialized by calling routine GRINIT. This routine asks about the type of 
terminal that Is being used. Second, routine CHRSIZ Is called to change the 
character size to #3 (next to the smallest) If a 4014- type terminal Is being 
used. Third, the screen Is erased, and the program asks the operator to 
supply four numbers which control the plot to be produced. The first value 
requested is N, which Is the nunber of points. N should range between 10 and 
30,000. The second value is the Shrinkage Factor (see below). This value 
should be between 0 and 2. The third number Is the Angle Increment, and the 
fourth Is the Line Type (normally 1, 97, 98, 99, 100, 104, or 112). The 
meaning of these parameters should be clear from the discussion following, and 
from observing the program's operation. The plot is drawn as follows. First, 
a point slightly to the right of the center of the terminal's screen Is the 
middle of the plot. From this point, any other point can be defined by a 
vector with a particular radius and angle, using polar coordinates. The 

program selects a starting value for the radius, and a starting angle of 0 
(relative to horizontal), so that the first point on the plot is to the right 
of the plot's center. The program then draws N lines, by changing the values 
for R and the angle each time. If the Shrinkage Factor Is 0, then R is always 
the same (all points will lie on an Imaginary circle). If the Shrinkage Is 1, 
then the last point will be at the center and all other points will He on an 
Imaginary spiral. Other values for the Shrinkage may be used also. After 
plotting each point on the plot, the program changes the angle by the user- 
specified Angle Increment. If this value Is very small, then a clrle or 
spiral will be drawn. If the value Is larger, then a polygon or star may be 
drawn. The process will continue until N lines have been drawn, and then the 
user may try a different combination of parameters to produce a new plot. 
Note that after R and the angle are computed for each point, they are 
transformed back Into rectangular (X and Y) coordinates before the call to 
MPLOT. The various line types are defined In the section on routine MPLOT. 
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Refer there for more Information. However, the various dotted and dashed 
lines (types 97-104) can only be produced on a Tektronix 4014 terminal with 
the enhanced graphics option. 

Although this program does not do anything particularly useful. It does 
demonstrate the use of several graphics subroutines. In particular. It Is an 
example of how MPLOT may be used In a special plotting situation, with little 
extra programming required. It also shows how a relatively simple program can 
generate rather Intricate plots. The program listing and sample output 
follow. 
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0001 

c 

c 

c 

n 

u 

IJ 

c 

c 

c 

c 

c 

c 

0002 

c 

c 

0003 

0004 

C 

0 

0005 

0006 

0007 

C 

0008 20 

0009 

0010 
0011 

C 

0012 

c 

0013 

0014 

C 

C 

C 

c 

0015 

0016 

0017 

0018 
('019 

0020 50 

0021 


PROGRAM DEMQGR 

DEMONSTRATION PROGRAM FOR GRAPHICS CAPABILITY 
OF TEKTRONIX TERMINALS 

C 0 MR I L I N G / L I NK I NG PROCEDURE ♦ 

FORT/WARN DEMOGR/LIST 

L INK/MAP :riKl/LIBlSY. ‘FPU/LIB JSYJWGCLIB DEMOGR 
NOTE: /LIBJSYJFPU SHOULD NOT BE USED WHEN YOU WILL 

BE RUNNING THIS PROGRAM ON A REGULAR LSI-11 
(USE IT ONLY FOR LSI-ll/23'S) 

RO = 1500. 

DEFINE CENTER OF DISPLAY PATTERN (SLIGHTLY TO THE RIGHT 
OF THE ACTUAL CENTER OF THE SCREEN) 

IXCNTR = 2500 
IYCNTR = 1550 

INITIALIZE GRAPHICS PARAMETERS 
CALL GRINIT(-1,0,0) 

CALL CHRSIZ<3) 

CALL ERASE 

TYPE *. 'Enter N ( 10-30000) r SHRINKAGE (0-2) >' 

TYPE 'ANGLE INCREMENT ( 1-360) » end LINE' 

TYPE * t ' TYPE (1»97»98»99»100»104»112) ' 

ACCEPT # »N> RF AC TR» ANGINC* I TYPE 

CHANGE ANGLE INCREMENT FROM DEGREES TO RADIANS 

ANGINC = ANGINC * 3.141593 / 180. 

MOVE TO RIGHT SIDE (STARTING POINT) 

CALL MPLOT ( IXCNTR+IFIX(RO) t IYCNTR»0 > 

DO 50 K=1 > N 

DRAW LINE TO A POINT AT RADIUS SLIGHTLY LESS 
THAN PREVIOUS VALUE (DETERMINED BY SHRINKAGE) 

AND AT AN ANGLE OF ANGINC COUNTERCLOCKWISE 
FROM THE PREVIOUS POINT 
ANGLE = FLOAT ( K ) * ANGINC 

R * RO * (1. - RFACTR * FLOAT ( K ) /FLOAT ( N ) ) 

IX = R * COS ( ANGLE ) + IXCNTR 
IY = R * S IN ( ANGLE ) + IYCNTR 
CALL MF'LOT ( IX» IY» I TYPE) 

CONTINUE 

CALL MPLOT (0»300f-l) 


0022 

0023 

0024 

0025 


: 'AUSE 'Hit 
CALL ERASE 

go ro 20 

END 


Return key' 
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ORIGilV L PASS V5 
OF POOR QUALITY 


FORTRAN IV 

Storage 

Map for Program Unit 

DEMOGR 



Local 

Variables* .PSECT 

♦DATA* 

Size 

= 000050 

( 20. 

words ) 


Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

ANGINC 

R*4 

000016 

ANGLE 

R*4 

000026 

I TYPE 

1*2 

000022 

IX 

I# 2 

000036 

IXCNTR 

I *2 

000004 

IY 

1*2 

000040 

IYCNTR 

I *2 

000006 

K 

1*2 

000024 

N 

1*2 

000010 

R 

R*4 

000032 

RFACTR 

R*4 

000012 

RO 

R*4 

000000 


Subroutines* Functions* Statement and Processor-Defined Functions* 


Name 

Type 

Name 

Type 

Name 

Type 

CHRSIZ 

R*4 

COS 

R'*4 

ERASE 

R#4 

IK IX 

1*2 

MPLOT 

1*2 

SIN 

R*4 


Name Type Name Type 

FLOAT R*4 6RINIT R*4 
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If. 


Enter N (10-30000). SHRINKAGE (0-2). 
ANCLE INCRERENT (1-360). m2 LINE 
TVPE (1.97.98.99.100.I04.U2) 

500, 0.5. 103,99 


?\\ S v VVnnN v\\ 




IBB 


»» 


Si; 


//AW 

w/A 


m 








"W 


»AOSE - Hit Return key '\wSTO 


The third sample program Is GRTEST. This demonstrates a more realistic 
plotting situation, although dunmy data was used for this example. The 
program uses GRINIT, CHRSIZ, GRID, ANOTAT, XYPLOT, MPLOT, COPY, and BELL. 
Several other routines are also called by these. Refer to the descriptions of 
each subroutine for more details. 

The operation of this program should be fairly clear from the program 
listing and the sample output. However, a few points are worth noting. 

First, vhen you use ANOTAT, you do not have to label every grid line. Every 
other vertical grid line was labelled In this example by specifying 10 
horizontal segments for GRID, and 5 for ANOTAT. Second, when drawing a single 
line (such as the "50% of max." line here), two calls to MPLOT are usually 
required. The first call, with the third argument set to 0, moves the current 
position to one end of the line (without drawing anything on the screen). The 
second call to MPLOT, with the third argument set to a positive number, draws 
the line to the coordinates of the specified end point. In addition, when 
labelling the plot axes or performing other alphanumeric I/O, you should first 
call MPLOT with the third argument set to -I, then use a formatted write (or 
type) statement, with a In the first print position In order to disable 

carriage control on that line. Otherwise, the line may be printed at a 

location different from the one you specified. 

Finally, the calls to routine BELL demonstrate how a program may provide 
some auditory feedback to the operator, perhaps to verify that certain data Is 
acceptable or unacceptable (using different sounds), or to let a user who Is 
away from the terminal i;now that some action Is needed. The program listing 
and sample output follow. 
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FORTRAN IV 


V02 .1-1 


PAGE 001 


OHKs. 

OF PC* ' ‘ 

W»d 23-Dec-81 15141146 


0001 


0002 


0003 


0004 


0005 


0004 


0007 


0008 


0009 


0010 


0011 

0012 

0013 

0015 

0016 

0017 


PROGRAM GRTEST 
C 

C DEMONSTRATE % TEST CERTAIN FEATURES OF TEKTRONIX 
C GRAPHICS PACKAGE 
C 

C COMPILING/LINKING PROCEDURE.* 

C FORT/WARN GRTEST/LIST 

C LINK /MAP J DK ♦ /LIB ♦ SY ♦ WGCLIB/LIB * SY ♦ FPU GRTEST 

C 

c note: the /lib:sy:fpu should not be used unless you 

C WILL BE RUNNING THE PROGRAM ON AN LSI-11/23. 

C DO NOT USE /LIB:SY:FPU WITH a REGULAR LSI-11. 

C 

c 

REAL X ( 30 ) , Y ( 30 ) 

c 

C DEFINE LIMITS OF PLOTTING AREA ON TERMINAL SCREEN 

DATA IL*IR,IB»IT / 400>3?00, 1000,2000 / 

C 

C DUMMY X & Y VALUES TO PLOT 

DATA X / 0. ,0.005,0.012,0.019,0.028,0.032,0.040,0.043, 

0 0 . 050 , 0 . 055 , 0 . 064 , 0 . 076 ,0.081,0. 095 ,0.100,0.104, 

0 0.109,0.116,0.122,0.125,0.134,0.143,0.151,0.158, 

0 0.166,0.170,0.175,0.183,0.190,0.195 / 

DATA Y / 0., 10. ,15. ,25. ,45. ,42. ,83. ,135. ,120. ,178. ,185., 

0 205. ,197. ,210., 225. ,222. ,265. ,308. ,332 .,322. ,345. , 

0 387. ,460. ,405. ,418. ,382. ,360. ,375. ,347. ,357. / 

C 

C INITIALIZE PARAMETERS FOR GRAPHICS ROUTINES (GET INFO. 

C FROM USER) 

CALL GRINIT (-1,0,0) 

C 

C SET CHARACTER SIZE TO #3 (NEXT TO SMALLEST) 

CALL CHRSIZ ( 3 ) 

C 

C DRAW HORIZ. & VERT. GRID LINES 

CALL GRID (10, 5, IL, I R, IB, IT, 97) 

C 

C LABEL (ANOTATE) GRID LINES WITH NUMERICAL USER UNITS 

CALL ANOTAT ( 5 , 5 , IL , IR , IB , I T ,0 . , 0 . 2 , 0 . , 500 . ) 

C 

C PLOT THE DATA (CONNECT POINTS WITH STRAIGHT LINES) 

CALL XYPLOT ( X , Y , 30 , IL , IR , IB , I T , 0 . , 0 . 2 , 0 , , 500 . , 1 > 

C 

C GET MAX. Y VALUE * DRAW HORIZ. LINE ON PLOT 

C AT LEVEL CORRESPONDING TO HALF THAT VALUE 

YMA X = 0. 

DO 50 K=1 , 30 

IF ( YMAX .LT. Y ( K ) ) YMAX = Y ( K ) 

50 CONTINUE 

IYHALF = ' t. 0 . 5* Y MAX ) / 500.; * (IT-IB.) + IB 
C MOVE TO LEFT SIDE OF PLOT AT PROPER Y DISTANCE UP 

CALL MPLPT ( IL. IYHALF, 0; 

C DRAW HORIZ. LINE TO RIGHT SIDE OF PLOT 
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0018 

C 

0019 

0020 

C 

0021 

002? 60 
C 

c 

c 

c 

c 

0023 
002 A 

0023 70 

C 
C 

0026 

C 

0027 

0028 

0029 80 
C 
C 

0030 

c 

c 

0031 

C 

c 

c 

0032 

0033 

0034 

C 

0035 

0036 


W*d 23-D»c-81 15141146 
CALL MPLOT < IR» I YHALF.99) 

MOVE TO POSITION JUST ABOVE THE LINE WE DREW 
IX « (IL+IR) / 2 + 250 
CALL MPLOT ( IX» I YHALF+10.-1 ) 

LABEL THE LINE 
TYPE 60 

FORMAT ('+50X of m«x.') 

NOTE THAT YOU MUST USE A + TO DISABLE CARRIAGE CONTROL 
IN THE FORMAT. OTHERWISE* A LINE FEED WOULD BE SENT TO 
THE TERMINAL BEFORE TYPING THE LINE. 

MOVE TO POSITION BELOW THE X-AXIS t LABEL IT 
CALL MPLOT ( IL+1350. IB-200. -1 ) 

TYPE 70 

FORMAT ('+STIMULUS DURATION Oumc.)") 

CHANGE TO CHARACTER SIZE *2 (SECOND LARGEST) 

CALL CHRSIZ (2) 

MOVE TO POSITION ABOVE TOP OF PLOT I LABEL IT 
CALL MPLOT < IL+600. IT+30. -1 > 

TYPE 80 

FORMAT < '+RESPQNSE AMPLITUDE VS. STIMULUS DURATION') 

MOVE TO BELOW BOTTOM OF PLOT 
CALL MPL0T(0. IB-300. -1) 

MAKE A HARD COPY OF THE TERMINAL SCREEN 
CALL COPY (1) 


MAKE SOME NOISES TO ALERT TJE OPERATOR IN CASE 
HE/SHE WENT TO SLEEP. 

CALI. BELL (20.8) 

CALL BELL (200.1) 

CALL BELL (5.60) 

STOP 

END 
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FORTRAN 10 Storase Map for Program Unit GRTEST 


Local 

Variables. .PSECT 

SDATA. 

Size 

- 000420 ( 

136. words) 


Name 

T VP9 

Offset 

Name 

Type 

Offset 

Name 

Twee 

Offset 

SB 

1*2 

000364 

IL 

1*2 

000360 

IR 

1*2 

000362 

IT 

1*2 

000366 

IX 

1*2 

000402 

IYHALF 

1*2 

000400 

K 

1*2 

000376 

YMAX 

R*4 

000372 





Local and COMMON Arrays 1 

Name Type Section Offset Size— Dimension* 

X R*4 SDATA 000000 000170 ( 60.) (30) 

Y R*4 $DATA 000170 000170 ( 60.) (30) 


Subroutines. Functions. Statement and Processor-Defined Functions. 


Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

ANOTAT 

R*4 

BELL 

R*4 

CHRSIZ 

R*4 

COPY 

R*4 

GRID 

R*4 

6RINIT 

R*4 

MPLOT 

1*2 

XYPLOT 

R*4 
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The fourth sample program is DISKRW. This routine demonstrates the use 
of subroutine DISKIO to read and write data to disk files. Refer to the 
program listing and sample output below. Note that this program Is written 
specifically for use with the DEC RT-11 operating system, and Is probably not 
adaptable to other systems. DISKIO calls a nunber of RT-11-speclfic subrou- 
tines, as Indicated In Its description, later In this doccment. DISKIO was 
written primarily to facilitate binary Input and output with variable record 
length to random (net necessarily sequential) blocks In a disk file, using a 
Fortran main program. The normal Fortran I/O normally requires that direct 
access binary (unformatted) records all be of the same length. Records In 

sequential files may be of variable length, but are Inconvenient and Ineffic- 
ient to access in a random-access fashion. 

In general, the normal DEC Fortran I/O methods, using an OPEN statement 
followed by READ or WRITE statements, are preferrable If their restrictions 
are not a problem. Otherwise, you may use subroutine DISKIO. However, 
records which have been written with DISKIO may only be read using the same 
routine, and not with the usual Fortran I/O methods. Basically, DISKIO allows 
more efficient I/O than Fortran in many cases, both In terms of time and also 
storage space, but its files are not compatible with those produced by the 
usual Fortran I/O routines. 

Operation of program DISKRW and DISKIO should be more clear after 
referring to the program listing, sample output, and output file dumps which 
follow. 


23 


r 
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0001 


VO 2 . 1-1 

PROGRAM DISKRW 


ORfcrAi r 

OF POOIK V . V 

Thu 01-0ct-81 15.12.58 


PAGE 001 


0002 

0003 

0004 

0005 

0006 
0007 


10 

20 

30 


PROGRAM TO VERIFY OPERATION i DEMONSTRATE USE OF SUBROUTINE BISKIO 
FOR WRITING DATA TO S READING IT FROM DISK FILES. 

SEE COMMENTS IN DISKIO FOR MORE INFO. 


WRITTEN BY} 
DATE? 


WILLIAM G. CROSIER 
11 JUNE 1980 


COMPILING/LINKING SEQUENCE} 

FORT/NOSWAP OISKRW/LIST 

FORT DISK I O/LI ST 

LINK/ MAP } DK } DISKRW. DISKIO 

INTEGER BUFR1 ( 1024 ) . BUFR2( 1024 ) . I ERR 
BYTE FILNAM< 12) 

DO 20 K= 1 » 12 
FILNAM(K)=0 
TYPE 30 

FORMAT ('OENTER COMPLETE FILE NAME FOR DISK I/O IN THE FORMAT'. 

8' 'DEVFILNAMEXT* . '/' WHERE} *DEV IS THE 3-CHARACTER DEVICE CODE'/ 
08X . ' * FILNAM * IS THE 6-CHARACTER BASIC FILE NAME'/ 
eaX.'^EXT* IS THE 3-CHARACTER EXTENSION/FILE TYPE DESIRED (OPT.)'/ 
0' EACH PORTION OF THE NAME SHOULD BE THE EXACT LENGTH SPECIFIED. V 
0' WITH SPACES ADDED. IF NECESSARY, TO FORM THE PROPER LENGTH.'/ 

OR PERIOD TO SEPARATE PORTIONS OF THE FILE', 
DK A12345DAT' / ' FILE NAME '? ',♦) 


0 ' 

@' 


DO NOT USE A COLON 
NAME.'/ ' EXAMPLE} 


f 


0003 


ACCEPT 40, FILNAM 

1 


0009 

40 

FORMAT (12A1) 

1 


0010 


TYPE ». 'ENTER SIZE OF FILE (IN NO. OF 256-WORD BLOCKS)' 

1 

- 

0011 


ACCEPT *, NBLK 

V 


0012 

45 

TYPE *, 'ENTER NO. OF WORDS TO WRITE S READ (1-1024)' 

i 

r 

7 ’ 

0013 


ACCEPT *, NWRDS 

5 

\ 

0014 


TYPE * » ' ENTER BLOCK NO. OF FILE WHERE I/O IS TO START' 



0015 


TYPE *,'(0=START AT 1ST BLOCK, 1=START AT 2ND, ETC.)' 


m ~ 

0016 


ACCEPT *, IBLK 



0017 


TYPE #, 'LEAVE FILE OPEN AFTER I/O (1=YES> ? ' 


* ‘ 

0018 


ACCEPT * , LOPEN 

* 


0019 


IWMODE = -1 


' * 

0020 


DO 50 K=l, NWRDS 

. 


0021 


BUFR1 (K)=K 

\ T 


0022 

50 

BUFR2 < K ) =0 

| 

• - 

■>023 


TYPE*, 'NOW WRITING DATA TO DISK' 


i 

0024 


CALL DISKIO (FILNAM, IWMODE, BUFR1 .NWRDS, IBLK, NBLK, IERR) 

r 


0025 


IF (IERR.NE.O) TYPE *, 'ERROR CODE' , IERR, ' DURING WRITE' 

} 



C FOR 

EXPLANATION OF ERROR CODES* SEE COMMENTS IN DISKIO 

1 

* ‘ 

0027 


IRMODE = 3 



0028 


IF (LOPF; .EQ. 1) IRMODE = -3 



0030 


TYPE * , ' NOW READING DATA FROM DISK' 


- 

0031 


CALL DISK 10 1 FILNAM , IRMODE , BUFR2 , NURDS , IBLK , NDUMM Y • IERR > 



0032 


IF (IERR.NE.O) TYPE *, 'ERROR CODE', IERR,' DURING READ' 



u()34 


TYPE 50, NWRDS • < BUFR2 ' K ) • K= 1 , NWRDS ) 



0035 

•50 

FORMAT ( ODATA READ FROM DISK ». SHOULD BE CONSECUTIVE INTEGERS' 


0/' FROM 1 THROUGH'. I! 


(1017) ) 
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Oh.w -A.. , 

OF POOR Vi’*-'"' • 4 


Thu 01-0ct-81 15512:58 

0036 IF (LOPEN . EQ . 1 ) GO TO 45 

0033 TYPE * « ' MORE (l=YES) ? ' 

0039 ACCEPT *. MORE 

0040 IF (MORE .EQ, 1) 60 TO 10 

0042 STOP 

0043 END 


FORTRAN IV Storage Map for Program Unit DISKRW 


local Variables* .PSECT $DATA» Size = 010050 ( 2068. words) 


Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

IBLK 

1*2 

010034 

I ERR 

1*2 

010024 

IRMODE 

1*2 

010042 

IWMODE 

1*2 

010040 

K 

1*2 

010026 

LOPEN 

1*2 

010036 

MORE 

1*2 

010046 

NBLK 

1*2 

010030 

N DUMMY 

1*2 

010044 

NURDS 

1*2 

010032 







Local 

and CQhMON Arrays 

♦ 

4 






N»me 

Type 

• Section 

Offset 

— 

— Size 

Dimensions 


BUFR1 

1*2 

♦ DATA 

000000 

004000 ( 1024.) 

(1024) 



FMJFR2 

1*2 

♦ DATA 

004000 

004000 ( 1024.) 

(1024) 



FILNAM 

L*1 

♦ DATA 

010000 

000014 ( 6.) 

( 12) 




Subroutines * Functions* Statement and Processor—Def ined Functions! 


Name 

T ype 

Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

DISK 10 

R*4 
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ORIGINAL F. r £-S 
Of POOR Q'L r u > 

DYl :« 00002 .DAT/N 
BLOCK NUMBER 00000 

000/ 000001 000002 000003 000004 000005 000006 000007 000010 

020/ 000011 000012 000013 000014 000015 000.016 000017 000020 

040/ 000021 000022 000023 000024 000025 000026 000027 000030 

060/ 000031 000032 000033 000034 000035 000036 000037 000040 

100/ 000041 000042 000043 000044 000045 000046 000047 000050 

120/ 000051 000052 000053 000054 000055 000056 000057 000060 

140/ 000000 000000 000000 000000 000000 000000 000000 000000 

160 / 000000 000000 000000 000000 000000 000000 000000 000000 

200/ 000000 000000 000000 000000 OCOOOO 000000 000000 000000 

220/ 000000 000000 000000 000000 000000 000000 000000 000000 

240/ 000000 000000 000000 000000 000000 000000 000000 000000 

260/ 000000 000000 000000 000000 000000 000000 000000 000000 

300/ 000000 000000 000000 000000 000000 000000 000000 000000 

320/ 000000 000000 000000 000000 000000 000000 000000 000000 

340/ 000000 000000 000000 000000 000000 000000 000000 000000 

360/ 000000 000000 000000 000000 000000 000000 000000 000000 

400/ 000000 000000 000000 000000 000000 000000 000000 000000 

420/ 000000 000000 000000 000000 000000 000000 000000 000000 

440/ 000000 000000 000000 000000 000000 000000 000000 000000 

460/ 000000 000000 000000 000000 000000 000000 000000 UOOOOO 

500/ 000000 000000 000000 000000 000000 000000 000000 000000 

520/ 000000 000000 000000 000000 000000 000000 000000 000000 

540/ OOCOOO 000000 000000 000000 000000 000000 000000 000000 

560/ 000000 000000 000000 000000 000000 000000 000000 000000 

600/ COOOOO 000000 000000 000000 000000 000000 000000 000000 

620/ 000000 000000 000000 000000 000000 000000 000000 000000 

640/ 000000 000000 000000 000000 000000 000000 000000 000000 

660/ COOOOO 000000 000000 OCOOOO 000000 000000 000000 000000 

700 / 000000 000000 000000 000000 000000 000000 000000 000000 

720/ 000000 000000 COOOOO 000000 000000 000000 000000 000000 

740/ 000000 000000 OOOCOO 000000 000000 000000 000000 000000 

760/ 000000 000000 000000 OOOOCO 000000 000000 000000 000000 
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SLOCK NUMBER 00001 

000/ 

000001 

000002 

020/ 

000011 

000012 

040/ 

000021 

000022 

060/ 

000031 

000032 

100/ 

OOOOOO 

oooooo 

170/ 

OOOOOO 

oooooo 

.140/ 

OOOOOO 

oooooo 

160/ 

oooooo 

oooooo 

200/ 

oooooo 

oooooo 

220/ 

oooooo 

oooooo 

240/ 

oooooo 

oooooo 

260/ 

oooooo 

oooooo 

300/ 

oooooo 

oooooo 

320/ 

oooooo 

oooooo 

340/ 

oooooo 

oooooo 

360/ 

oooooo 

oooooo 

400/ 

oooooo 

oooooo 

420/ 

oooooo 

oooooo 

440/ 

oooooo 

oooooo 

460/ 

oooooo 

oooooo 

500/ 

oooooo 

oooooo 

520/ 

oooooo 

oooooo 

540/ 

oooooo 

oooooo 

560/ 

oooooo 

oooooo 

600/ 

oooooo 

oooooo 

620/ 

oooooo 

oooooo 

640/ 

oooooo 

oooooo 

660/ 

oooooo 

oooooo 

700/ 

oooooo 

oooooo 

720/ 

oooooo 

oooooo 

740/ 

oooooo 

oooooo 

760/ 

oooooo 

oooooo 


000003 000004 000005 
000013 000014 000015 
000023 000024 000025 
000033 000034 000035 

000000 000000 oooooo 
000000 oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 
oooooo oooooo oooooo 


000006 000007 
000016 000017 
000026 000027 
000036 000037 

oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 
oooooo oooooo 


000010 

000020 

000030 

000040 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 

oooooo 
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GRAPHICS 

SUBROUTINES 
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Subroutine GRINIT 


This routine Initializes the various parameters used by other graphics 
routines, and should be called before any of the other graphics routines. It 
Is not needed If only the data acquisition routines ISAMPA, DAQISR, DAQ or 
miscellaneous routines BELL, WAIT, or DISKIO are used. 

Information on the parameters or argunents to be passed to GRINIT appears 
In the program listing below. When writing a particular program. If you will 
always be using It with the same type of terminal and hard copy unit, then you 
may specify their model nunbers In the argument list so that GRINIT may define 
Its Internal parameters. However, If you may use different types of 
terminals, then you should set the first argument to -1 so that GRINIT will 
query the user at execution time about the type of terminal and hard copy unit 
to be used. The program listing follows. 
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1 


ORICU’J 


OF POOR Q !5 ^ '‘ ' 


FORTRAN IV 
0001 


0002 

0003 


0004 

0005 

0006 
0007 


000 ? 

0010 


V02.1-1 Wed 23-Oec-Sl 10:2lM8 
SUBROUTINE GRINIT ( ITERM* IHDCP* IENH) 


PAGE 001 


0014 

0015 

0017 

0018 


001 ? 

0020 

0021 

0022 


■>024 


WILLI Art G. CROSIER 
28 SEPT. 1981 


L WRITTEN BY: 

c revised: 
c 

il purpose: initialize parameters defining terminal 

C CHARACTERISTICS FOR THE VARIOUS GRAPHICS ROUTINES 

c arguments: 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


I TERM = MODEL NUMBER OF TERMINAL BEING USED. 

IF =-l* THEN USER WILL BE QUERIED DURING PROGRAM 
EXECUTION ABOUT TERMINAL i HARD COPY UNIT. 

IHDCP = MODEL NUMBER OF HARD COPY UNIT BEING USED 
SET IHDCP=0 IF NO HARD COPY UNIT IS BEING USED 
IENH = FLAG FOR ENHANCED GRAPHICS OPTION FOR TERMINAL. 
SET IENH=1 IF YOU HAVE A 4014/4015 WITH THIS OPTION. 
OTHERWISE* SET IENH=0 


NOTE! IHDCP S IENH ARE BOTH IGNORED IF ITERM=-1. IN THIS CASE* 
THEN THE USER WILL BE ASKED TO SUPPLY THE APPROPRIATE 
INFO. FOR THE TERMINAL WHEN THIS ROUTINE IS CALLED. 


0 


L0GICAL*1 ANSWER 

COMMON /GRFCOM/ MCHRSZ*LSIZE* I WIDTH* IHIGHT * IENHAN* 
TIMERA* riMHDC 


I TERM 1 = ITERM 
IHDCF'l = IHDCP 
IENH1 = IENH 
IF (ITERM1 .GT. 0) 


60 TO 60 


20 


0011 

0012 

0013 


30 


40 


QUERY USER ABOUT TYPE OF TERMINAL % HARD COPY UNIT BEING USED 
TYPE 20 

FORMAT (' What is the model no. for the terminal sou are '* 

0 'using? <Don''t enter' / ' dash numbers. Example 

(? ' If you have a 4014-1* Just enter 4014) ?'*$) 

ACCEPT ** ITERM 1 
TYPE 30 

FORMAT (' What is the model no. for the hard copy unit'* 

0 ' you are using?' / ' (If none is connected to'* 

0 ' your terminal* enter 0.) ?'•$) 

ACCEPT #* IHDCF’l 

IF aTERMl.NE.4014 .AND. I TERM1 . NE . 4015 ) 60 TO 60 

TYPE 40 

FORMAT (' Does your terminal have the enhanced graphics'* 

0 ' option?' / ' (Can it drau dotted and dashed'* 

0 ' lines?) (Y=Yes; '**> 

ACCEPT 50* ANSWER 
FORMAT (Al> 

IENH1 = 0 

IF i ANSWER . £Q . ' Y ' ) IENH1 = 1 


C 

60 


MCHRSZ = 1 

IF USING A 4014/4015 TERMINAL* MULTIPLE 


IARACTEF: SIZES ARE AVAIL. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


FORTRAN 10 002.1-1 Wed 23-Dec-81 10J21J48 PAGE 002 

0025 IF < ITERM1 .NE.4014 .AND. ITERM1 . NE . 4015 > fICHRSZ * 0 

C 

C GET TINE TO ALLOW FOR ERASING SCREEN * 1,5 SEC, 


0027 

TIMERA = 1.5 




C SET DEFAULT HARD CCPY TIME - 20 SEC. 



0029 

TIMHDC = 20.0 



0029 

IF ( IHDCP1 ,LE. 0) TIMHDC=0 

! NO 

HARD COPY UNIT AVAIL. 

0031 

IF ( IHDCP1 .EG. 4631) TIMHDC-10.0 

! 10 

SEC FOR MODEL 4631 

0033 

IENHAN = IENH1 



0034 

RETURN 



0035 

END 




C 


FORTRAN 10 Storage Nap for Program Unit GRINIT 

Local Variables. .PSECT *DATA» Size * 000016 ( 7. words) 

Name Type Offset Name Type Offset Name Type Offset 

ANSWER L*1 000006 IENH 1*2 9 000004 IENH1 1*2 000014 

IHDCP 1*2 9 000002 IHDCP1 1*2 000012 ITERN 1*2 S 000000 

ITERN1 1*2 000010 


COMMON Block /GRFCOM/ » Size = 000022 ( 9. words) 

Name Type Offset Name Type Offset Name Type Offset 
MCHRSZ 1*2 000000 LSIZE 1*2 000002 IWIDTH 1*2 000004 
IHIGHT 1*2 000006 IENHAN 1*2 000010 TIMERA R*4 000012 
TIMHDC R*4 000016 
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Subroutine MPLOT 


This Is the basic line-drawing subroutine. It Is used by all of the 
other graphics routines which draw lines or move the cursor, and It may be 
called directly by the user's program also. 

The argunents or parameters for this routine are described In the program 
listing below. Refer there for more Information. In addition, there are 
further details, concerning the various line types available In a 4014 with 
enhanced graphics, In the Tektronix manuals for the 4014 terminal. If you are 
not using a 4014 or 4015, then GRINIT will cause all requests for line types 
other than 1 (normal solid line) to be responded to as If line type 1 was 

specified. In this way, you do not need to worry about the presence of the 

optional enhanced graphics hardware In your terminal when you write the 
software, because If the hardware Is not present then solid lines will be 
drawn regardless of the line type you specify with parameter IPEN. 

Besides drawing a line, MPLOT can also be used to move the cursor 

Immediately before plotting by setting parameter IPEN to 0. Also, It can be 

used to move the cursor before typing alphanumeric data by setting IPEN to -1. 
Refer to the third sample main program (GRTEST) In an earlier section of this 

report for examples of this. The second sample pgoram (DEMOGR) may also be 
helpful when using MPLOT. 

If you are not using a 4014/4015 terminal with enhanced graphics, then 
you will have only 1024 point horizontal resolution horizontally and 780 point 
resolution vertically, but the same coordinates will still refer to the same 
point on the screen. Thus, an X coordinate of 2048 and a Y coordinate of 1560 
will always refer to the center of the screen, regardless of the type of 

Tektronix terminal you are using. 

If you want to modify this routine for use with another type of terminal 
or with an X-Y plotter, then several parts of the subroutine will need to be 
changed. First, Tektronix requires that the X and Y coordinates be split up 
and sent In a particular sequence. This sequence, under the heading "Draw 
Vector" In the program listing, would probably need to be changed for a 
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dlffe-ent terminal or plotter. Second, many plotters require a delay after 
transmitting the coordinates for a new point so that the pen has time to move. 
For this purpose, you can use a call to subroutine WAIT. If an analog X-Y 
recorder was to be used rather than a digital plotter, then the X and Y 
coordinates would need to be sent to the plotter via two IPOKE calls (refer to 
the System Subroutine Library routines) to a pair of digital -to-analog 
converters. The outputs of these may have to be slowed down with a pair of 
matched R-C networks In order to avoid too^rapld changes In the pen position, 
If the recorder so requires. Using this routine with a digital plotter should 
be much simpler, however. Many digital plotters also have built-in character 
generators for drawing the standard ASCII character set also, while analog X-Y 
recorders do not have this capability. 

Note that when using a line type of 112 (for the write-thru mode), the 
lines which are drawn are faint, and do not store on the CRT screen. Norm- 
ally, when using this line type you should refresh the display by repeating 
the plotting of the appropriate lines. Ideally, the lines should be redrawn 
at least 50 times per second In order to prevent flickering and to make the 
display easier to see. However, slower refresh »*ates may be necessary If a 
large mmber of lines have to be redrawn with each repetition. In addition, 
static displays (with line types other than 112) may be combined with dynamic 
displays (line Type 112) at the same time. That way you only have to refresh 
the lines whose positions change. 

One problem when using MPLOT (and any routines which call It) has 
appeared In some systems when using the RT-11 Version 4 Foreground/Background 
monitor. The problem appears to be In the RT-11 terminal handler (TT.SYS) or 
In the System Subroutine Library routine ITTOUR, since the difficulty has 
never appeared with the RT-11 Single Job monitor. If your plots do not come 
our correctly but appear garbled, try running your program under the Single 
Job monitor (RT11SJ) rather than the Foreground/Background one (RT11FB). 

MPLOT requires the following routine from the System Subroutine Library: 
ITTOUR 

A listing of MPLOT follows. 
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FORTRAN IV 


ORIGINAL PAC-S IS 

OF Por-> r ' vy 

W*d 23-D«c-8l 10:21:59 


PAGE 001 


0001 


SUBROUTINE MPL0T( IXf IYf IPEN) 


C- 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


author: william g. crosier 

REVISED 25 SEPT. 1981 

BASED ON ROUTINE WRITTEN BY CHUCK MANN 


TECHNOLOGY INCORPORATED 


IPEN 

IPEN 


PURPOSE: TO DRAW LIGHT OR DARK VECTORS FROM THE CURSOR'S 
(PEN'S) PRESENT LOCATION TO THE COORDINATES PASSED. 

arguments: 

IX » X COORDINATE IN TEKTRONIX SCREEN UNITS (0 TO 4095) 
IY * Y COORDINATE IN TEKTRONIX SCREEN UNITS <0 TO 3120) 
IPEN CONTROLS PLOTTING AS FOLLOWS: 

IPEN = 0 MOVES POSITION TO aXtlY) WITHOUT 
DRAWING A LINE ( DAiAK VECTOR) 

0 MOVES POSITION * CHANGES TO ALPHA MODE 
0 DRAWS A VISIBLE LINE TO <IX»IY) 

IF IPEN * If THEN NORMAL SOLID LINE 
IF IPEN * 97 f THEN DOTTED LINE (ON 4014) 

» 98 f THEN DOT-DASH LINE 

* 99 f THEN SHORT DASH LINE 
» 100f THEN LONG DASH LINE 

* 104 r THEN DEFOCUSED Z-AXIS 
(SLIGHTLY WIDER LINE) 

■ 1 12f THEN WRITE-THRU (NON-STORE) 

(THE VARIOUS DOTTED t DASHED LINES 
WILL ONLY BE PRODUCED ON A 4014 



c 

c 


TERMINAL WITH ENHANCED GRAPHICS 

OPTION. ) 

0002 

c 

COMMON /GRFCOM/ MCHRSZf LSIZEf IWIDfHf IHlGHTf 

IENHANf 


r* 

8 

TIMERA 

f TIMHDC 


0003 

L 

DATA LSIZEf 

IWIDIH* 

IHIGHTf IENHANf TIMERAf TIMHDC 



& 

/ If 56 f 88 f Of 10. Of 1 .5 / 


0004 


DATA MINX.MAXXfMINYfMAXY /O r 4095 f 0 f 3 120/ 


0005 


DATA IPEN1 

/O/ 



0006 


1X1 = IX 




0007 


I r 1 = IY 




0008 


IF (IX .LT. 

MINX) 

IX1=MINX 


0010 


IF (IX .GT. 

MAXX ) 

IX1-MAXX 


0012 


IF (IY .LT . 

MINY > 

IY1*MINY 


0014 


IF (If .GT. 

rtAXY ) 

I Y 1 *MAX Y 


0016 

c 

c 

IF ( IPEN.GT 

. 0 > GO 

TO 10 



DRAW DARK VECTOR (MOVE POSITION) 


00 1 S 

100 

IF (ITTOUR( 

*35) .NE 

.0) GO T 0 100 

! SEND GS 

0020 

10 

IF (IPEN.LE 

* 0 * OR 4 

IENHAN.EQ.0 .OR. IPEN. EG 

. I PEN 1 > 


c 

£ 

GO TO 

30 



35 


t 


ft. 


FORTRAN IV 


C 

C 


0022 

0023 

0025 

20 

0027 

25 

0029 

30 

0030 

C 


C 


c 

0031 

0032 

Lr 

105 


c 

0034 

0035 

108 


c 

0037 

0038 

110 


c 

0040 

0041 

115 


C 

0043 

0044 

120 

0046 

C 

r 

0048 

1 v 

200 

0050 

900 

0051 

c 


ORIGINAL PAGE \3 
OF POOR O' ' A' !TY 


V02.1-1 Wed 23-£i»c-81 10J21 J59 


PAGE 002 


CHANGE TO SELECTED PLOTTING MODE IF USING TERMINAL WITH 
ENHANCED GRAPHICS 


LINTYP * IPEN 

IF (IPEN .LT. 96) LINTYP * 96 
IF < ITT0UR<27) .NE.O) GO TO 20 
IF (IT TOUR (LINTYP) .NE.O) GO TG 25 
CGNTINUE 
IPEN1 » IPEN 


! SOLID LINE? 
! SEND ESC 
! SET LINE TYPE 

(RESET PREV ♦ 


VALUE 


DRAW VECTOR 


SEND HIGH ORDER Y BYTE 
ICH**40+IY1/128 

IF (ITTOUR(ICH).NE.O) GO TO 105 
SEND EXTRA BYTE (2 LSB'S OF X t Y ) 

ICH* * 140 + ( (IY1.AND. *3>*4) + <1X1. AND. *3) 
IF < ITTQUR< ICH) .NE. 0) GO TO 108 
SEND LOW ORDER Y BYTE 
ICH* * 140+ < (IY1/4) .AND. *37) 

IF (ITT OUR (ICH). NE.O) GO TO 110 
SEND HIGH ORDER X BYTE 
ICH* *40+ 1 XI/ 128 

IF ( I T TOUR ( ICH ) , NE . 0 ) GO TO H5 
SEND LOW ORDER X BYTE 
ICH»*100+< (1X1/4) .AND. *37) 

IF(ITTOURdCH) .NE.O) GO TO 120 
IF (IPEN .GE. 0) GO TO 900 


CHANGE TO ALPHA MODE IF IPEN IS LESS THAN 0 

IF (ITT0URC37) .NE. 0) GO TO 200 

RETURN 

END 


FORTRAN IV Storage Map for Program Unit MPLOT 


Local 

Variables* .PSECT 

♦ DAT A f 

Size 

* 000030 ( 

12. 

words ) 


Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

ICH 

1*2 

000026 

IPEN 

1*2 9 

000004 

IPEN1 

1*2 

000016 

IX 

1*2 8 

000000 

1X1 

1*2 

000020 

IY 

1*2 8 

000002 

IY1 

1*2 

000022 

LINTYP 

1*2 

000024 

MAXX 

1*2 

000010 

MAXY 

1*2 

000014 

MINX 

1*2 

000006 

MINY 

1*2 

000012 


COMMON 

Block 

/GRFCOM/ » 

Size = 

000022 

( 9. 

words > 



Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

MCHRSZ 

1*2 

000000 

LSIZE 

1*2 

000002 

IUIDTH 

1*2 

000004 

IHIGHT 

TIMHDC 

1*2 

R*4 

000006 

000016 

IENHAM 

1*2 

000010 

TIMERA 

R*4 

000012 

Subroutines 9 

Functions * 

Statement and 

Processor- Defined 

Functions : 


Name Type Name Type Name Type Name Type Name Type 
ITTOUR 1*2 
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Subroutine COPY 


I 


This short routine causes a hard copy unit (If one is present) to make a 
p*per copy of whatever appears on the Tektronix terminal screen. After 
waiting for a fixed length cf time, or until the user strikes the Return key 
(depending on the parameter IFLAG), the routine returns to the calling 
program. More details are given In the program listing. 

COPY requires the following subroutines: 

ERASE 

WAIT 

ITTOUR (from the System Subroutine Library) 

A listing of routine COPY follows. 
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ORIGINAL L. - h ' 

OF POOR QUALIlY 


l- OR TRAN IV 


V02.1-1 Wed 23-Dec-P 1UJ22J12 


PAGE 001 


0001 


0002 


0003 

0005 

0007 


0009 

0010 
0012 


0014 


0015 

0016 
0017 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE COPYaFLAG) 


PURPOSE: TO HAKE a HARD COPY OF THE TEKTRONIX SCREEN. 

if: iflag=o» waif for user to type <cr>. then 

ERASE SCREEN. 

= 1* RETURN AFTER MAKING COPY 1 DON ' T 
WAIT FOR RESPONSE r DON'T ERASE) 

DON'T WAIT FOR RESPONSE AFTER MAKING 
COPY t BUT ERASE SCREEN. 


C 

C 

COMMON /GRFCQM/ MCHRSZ . LSIZE . I WIDTH • IHIGHT » IENHAN t 
0 TIMERAf TIMHDC 

C 

C MAKE HARD COPY OF THE SCREEN IF HARD COPY UNIT IS AVAIL. 

IF (TIMHDC .LE. 0) GO TO 50 

230 IF (ITT OUR ( 27 ) . N£ . 0 ,■ GO TO 230 (SEND ESC 

235 IF ( ITT0UR(23) .NE.O) GO TO 275 (SEND ETB 

C 

C WAIT FOR SCREEN TO COPY 

CALL WAIT ( TIMHDC » 0 ) (WAIT TIMHDC SEC. 

50 IF (IFLA6.EQ.1) RETURN 

IF ( IFLAG . EQ . 2 ) GO TO 200 
C 

C WAIT FOR USER TO RESPOND 

PAUSE. 'HIT RETURN KEY TO CONTINUE' 

C 

C CLEAR THE SCREEN 
200 CALL ERASE 
RETURN 
END 

C 


FORTRAN IV Storsse Mae for Program Unit COPY 


Locs 1 

Variables. .PSECT SDATA. 

Size 

= 000002 

v 1* words) 


Name 

Type 

Offset 

Name 

T y^e 

Offset 

Name 

Type 

Offset 

IFLAG 

1*2 i ? 

000000 







COMMON 

Block 

/ GRFCOM/ . 

Size = 

000022 9. 

words ) 



Name 

T ype 

0 f f s e t 

Name 

T ype 

Offset 

Name 

T ype 

Of f set 

mCHRSZ 

1*2 

000000 

LSIZE 

1*2 

000002 

IWIDTH 

1*2 

000004 

IHIGHT 

1*2 

000006 

IENHAN 

1*2 

000010 

TlhERA 

K*4 

0000 1 2 

TIMHDC 

R* 4 

OOOOlo 







^ ijbroux i nes * 

Functions 

t Statement a 

no Processor-Defined 

Functions * 

Nc me 

T ype 

Name T 

ype Ncir.e 

Tv?e Ns 

me Typ*? 

Name 

T y p e 

ERASE 

R'*4 

IT TOUR 

1*2 WAIT 

R*4 





Subroutine ERASE 


This simple routine Is used to erase or clear the Tektronix terminal 
screen. It sends an escape followed by a form feed character, waits a short 
time for the erase process to finish, then returns to the calling program. 

ERASE requires the following subroutines: 

WAIT 

ITTOUR (from the System Subroutine Library) 

A listing of routine ERASE follows. 
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ORIGINAL PAC2 m 
OF POOR DUALITY 

FORTRAN IV V02.1-1 Wed 23-Dec-81 10:22121 PAGE 001 

0001 SUBROUTINE ERASE 
C 

C ERASE (CLEAR) THE SCREEN ON THE TEKTRONIX TERMINAL 

C 

0002 COMMON /GRFCOM/ HCHRSZfLSIZEr I WIDTH* IHIGHT * IENHAN* 

<? TIMERA * TIMHDC 

C 

0003 10 IF ( I T TOUR ( 27 ) .NE. 0) GO TO 10 ! SEND ESC 

0005 20 IF ( I T TOUR ( 12 ) .NE. 0) GO TO 20 ! SEND FF 

0007 CALL WAIT(TIMERAfO) ! WAIT TIMERA SEC. 

C 

0008 RETURN 

0009 END 
C 



FORTRAN 10 Storage Map for Program Unit ERASE 

COMMON Block /GRFCOM/* Size = 000022 < 9. words) 

Name Type Offset Name Type Offset Name Type Offset 

MCHRSZ 1*2 000000 LSIZE 1*2 000002 IWIDTH 1*2 000004 

IHIGHT 1*2 000006 IENHAN 1*2 000010 TIMERA R*4 000012 

TIMHDC R*4 000016 

Subroutines* Functions* Statement and Processor-Defined Functions: 

Name Type Name Type Name Type Name Type Name Type 

ITTOUR 1*2 WAIT R*4 
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Subroutine CHRSIZ 


This subroutine can be used to set the character size for those Tektronix 
terminals with multiple sizes available (the 4014/4015). This routine should 
be called before the other graphics routines which print alphanumeric 
characters, such as ANOTAT (or XYPLOT with ICOOE equal to 0). It can also be 
used before regular type statements In order to change to smaller or larger 
characters for page headings, etc. Once CHRSIZ Is called, the new character 
size stays in effect until CHRSIZ Is called again with a different size 
specified, even If the screen Is erased or If you change between alpha and 
graph modes. The parameter/ argument ISIZE sets the character size as 
described In the program listing. 

CHRSIZ requires the following subroutines from the System Subroutine 
Library: 


ITTC'JR 

A listing of CHRSIZ follows. 
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FORTRAN 10 


002 . 1-1 


PAGE 001 


OF POOR 


^ r- r' tf *7 


QUA 


Wed 23-Dec-81 10 : 22 J 30 


0001 

C 

C 

C 

c 

c 

c 

c 

c 

0002 

0003 

0004 


0005 

0006 
0007 

c 

0003 
0010 20 

0012 30 

0014 

0015 

0016 

C 

0017 99 

0013 

C 


SUBROUTINE CHRSIZ( ISIZE) 

ROUTINE TO SET CHARACTER SIZE ON TEKTRONIX 4014 SCREEN 

ARGUMENT ISIZE CONTROLS CHAR. SIZE AS FOLLOWS: 

ISIZE = l SELECTS LARGEST CHARACTERS 
ISIZE = 2 SELECTS HEDIUH-LARGE CHAR. 

ISIZE = 3 SELECTS MEDIUM-SMALL CHAR. 

ISIZE = 4 SELECTS SMALLEST CHAR. 

L0GICAL*1 ICODE ( 4 ) 

INTEGER IS IZE* IW< 4 ) * IH< 4 ) 

COMMON /GRFCOM/ MCHRSZ * LSIZE * IUIDTH * IHIGHT * IENHAN * 

<? TIMERA * TIMHDC 

DATA MCHRSZ / 1 / 

DATA ICODE / '8' * '9' * ' J ' * ' * ' / 

DATA IW / 56*51*34*31 /*IH / 88*82*53*48 / 

IF (ISIZE. LT.l .OR. ISIZE. GT. 4 .OR. MCHRSZ. EG. 0) GO TO 99 
IF < ITTOUR ( 27 ) .NE. 0) GO TO 20 ! SEND ESCAPE 

IF ( ITTOUR( ICODE(ISIZE) ) .NE. 0) GO TO 30 ! SEND CHAR. 

IUIDTH = IU( ISIZE) 'WIDTH OF A CHAR. 

IHIGHT = IH< ISIZE) (HEIGHT OF A CHAR. 

LSIZE = ISIZE 

RETURN 

END 


FORTRAN 10 Storage Map for Frog ram Unit CHRSIZ 


Local Oariables * .PSECT iDATA* Size = 000026 ( 11. words) 


Name 

Type 

Offset 

Name 

Type Offset 

Name 

Type Offset 

ISIZE 

1*2 £ 

000000 






COMMON 

E'lock 

/ GRFCOM/ * 

Size = 000022 

( 9. 

words ) 


Name 

Type 

Offset 

Name Type 

Offset 


Name Type 

Offset 

MCHRSZ 

1*2 

000000 

LSIZE 1*2 

000002 


IUIDTH 1*2 

000004 

IHIGHT 

1*2 

000006 

IENHAN 1*2 

000010 


TIMERA R*4 

000012 

TIMHDC 

R*4 

000016 






! ocel 

arid COMMON Arrays 

♦ 

* 





Name 

Typ 

e Section 

Offset 

— Size — 

— 

Dimensions 


ICODE 

L*1 

JDA7A 

000002 000004 ( 

2* ) 

< 4 > 


I’M 

1*2 

iDATA 

0000 lo 000010 ( 

4 ♦ > 

< 4 > 


IW 

1*2 

*DATA 

000006 000010 ( 

4* ■ 

( 4 ,» 



Subroutines* Functions* Statement and Processor-Defined Functions: 

Name Type Name Type Name Type ns me Type Name Type 

ITTOUR 1*2 
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Subroutine GINPUT 


This routine Is used with those terminals such as the 4010 and 4014, 
having user-controlled crosshairs for Interactive graphics Input fro* the 
user. Each time this routine Is called, the terminal switches to Graphics 
Input mode and displays the user-controlled crosshairs. The position of these 
two lines (one horizontal and one vertical) may be changed by the operator by 
turning a pair of thumbwheels by the keyboard. The crosshair lines are dim 
and do not store on the terminal screen. By using this routine, the program 
can allow the user to select a particular point on the display (perhaps one 
previously plotted with ARYPLT or XYPLOT) and can then easily get the 
coordinates of that point from the user. The subroutine does this by waiting 
until the user has positioned the crosshairs as desired. Then he or she can 
strike a single key on the keyboard and then the Return key. When this Is 
done, the terminal automatically transmits the first character that was struck 
along with the X and Y coordinates of the crosshairs to the computer. 
Subroutine GINPUT then returns the ASCII equivalent of the character and the 
crosshair coordinates to the calling program through the argument list. 

Each time GINPUT Is called, a single character and the coordinates of one 
point are sent to the program. The character struck does not appear on the 
screen. (The terminal "bypass" circuitry keeps this from happening.) The 
character transmitted can be used as a code for the program to Instruct It 
what function should be performed next, or It can be used for any other 
purpose by the calling program. Although both the X and Y coordinates are 
always passed to the program, frequently only one or the other Is needed. In 
that case the user can be Instructed to Ignore the crosshair line that Is not 
needed (for example, the horizontal one), and use only the other one to pick 
out the feature of Interest on the screen. This can speed up the program's 
operation by cutting In half the amount of time the operator must spend In 
positioning the crosshairs. This Is especially useful when one Is selecting a 
number of points from a plot of sampled analog data (such as may be produced 
by ARYPLT) versus time, for Instance, since only the vertical crosshair need 
be used to point out the times of Interest. The program can then look up the 
Y coordinate from the sampled data point In memory which corresponds to that 
time or X coordinate. 
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Although the graphical output routines In this package can plot points 
with 4096 point resolution If a 4014/4015 terminal with the enhanced graphics 
option Is used, the Tektronix terminals are all limited to 1024 point 
resolution for graphical Input . Therefore, In detailed plots, It may be 
Impossible to always select the exact point of Interest with the GINPUT 
routine. In general, your program should take the coordinates returned by 
GINPUT and find the point which you plotted that Is closest to the returned 

coordinates. Even then. It may be Impossible to always select a single given 
point from among several very closely spaced points. Usually, however, the 
error In selecting among closely spaced points Is not noticeable. 

Routine GINPUT requires the following subroutine from the System 
Subroutine Library: 

ITTOUR 

The listing for GINPUT follows. 
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1 


ORIGINAL FACE IS 
OF POOH QUALITY 


FORTRAN 10 


002 . 1-1 


Wed 23-Dec-Ql 10J22J4? 


RAGE 001 


OOOi 


0002 


0003 

0005 


L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

20 

30 

C 

c 

c 

c 


SUBROUTINE G INPUT ( I CHAR » IX » IY ) 

SWITCH TO GRAPHIC INPUT MODE* DISPLAY USER-CONTROLLED 
CROSSHAIRS. WAIT FOR USER TO TYPE IN A CHARACTER* 
t RETURN THAT CHARACTER ( I CHAR ) * AND THE X AND Y 
COORDINATES OF THE CROSSHAIRS TO THE CALLING 
PROGRAM . 

INTEGER I CHAR .HIGHX *LOUX*HIGHY*LOWY 

SWITCH TO GRAPHIC INPUT MODE * DISPLAY CROSSHAIR 


IF ( I TTOUR ( 27 ) .NE. 0) GO TO 20 
IF ( ITTOUR ( 26 ) .NE. 0) GO TO 30 


(SEND 
! SEND 


ESCAPE 

SUB 


0007 

ACCEPT 

0008 50 

FORMAT 

0009 

IX = 4 

0010 

IY = 4 

0011 

RETURN 

0012 

END 


WAIT UNTIL USER HITS A TERMINAL KEY * THEN GET THAT 
CHARACTER AND THE CROSSHAIR ADDRESS 


50 * I CHAR * HIGHX * LOWX » HIGHY * LOWY 


* < ( (HIGHX. AND. *37) * 

* < < (HIGHY.AND. *37) * 


32) + (LOUX.AND. *37) 
32) + (LOWY. AND. *37) 


FORTRAN IV Storaae Map for Program Unit GINF'UT 


Local Variables* .F'SECT SDATA* Size = 000016 ( 7. words) 


Name 

Type 

Offset 

Name 

Type Offset 

Name 

Type 

Off set 

HIGHX 

1*2 

000006 

HIGHY 

1*2 000012 

I CHAR 

1*2 0 

000000 

IX 

1*2 (? 

000002 

IY 

1*2 0 000004 

LOWX 

1*2 

000010 

LOWY 

1*2 

000014 






Subroutines 9 

Functions* Statement and Processor 

-Def ined 

Functions • 

Name 

Type 

Name 

Type N 

ame Type N3me 

T ype 

Name 

Type 

ITTOUR 

1*2 
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Subroutine ARYPLT 


This routine can be used to plot an Integer array of Y values with equal 
Increments In the X (horizontal) direction. It Is particularly well suited 
for plotting sampled data such as from an analog-to-dlgltal converter, with 
the samples equally spaced In time. All points are connected with straight 
solid lines. For Integer data, ARYPLT provides a significant savings In 
memory usage over XYPLOT, at least If a large nuaber of points are plotted. 
This Is so because real arrays In PDP-11 Fortran require two words to store 
each element, while Integer arrays only require one word per element. In 
addition, a separate array of X coordinates must be passed to XYPLOT, but Is 
not needed for ARYPLT since ARYPLT generates the X coordinates automatically. 

The arguaents or parameters are described briefly In the program listing. 
In addition, you should note that all of the arguments except for YSCALE are 
of Integer t>pe. If IYOFST Is set to 1560, then positive values In IARRAY 
will be plotted above the middle of the screen and negative values will be 
plotted below the middle. IYOFST can be set larger or smaller In crder to 
move the plot up or down, respectively. YSCALE can be set to a value of 1.0 
for many applications, If the elements of IARRAY are within a range of about 
-1000 to +1000, but YSCALE can be changed to adjust the vertical scale factor 
of the plot. 

Several calls to ARYPLT can be made without erasing the screen, with 
different values for IYOFST and/or LEFT and RIGHT, In order to plot several 
curves on the screen together. 

Routine ARYPLT requires the following subroutine. In addition to those 
from the Fortran library: 

MPLOT 

A listing of ARYPLT follows. 
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ORIGIN WL PAGE IS 
OF POOR QUALITY 

FORTRAN 10 002.1-1 Wed 23-Dec-81 10:22:53 PAGE 001 


0001 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0002 

0003 

0004 

0005 

0006 

c 

0007 

c 

0008 

0009 

0010 

0011 

0012 20 

0013 

0014 


SUBROUTINE ARYF'LT ( I ARRAY ?N? IYOFST ? YSCALE ? LEFT ? RIGHT ) 

ROUTINE TO PLOT AN INTEGER ARRAY OF Y VALUES 
IN * I ARRAY * ON THE TEKTRONIX TERMINAL. 

THE X-COORDINATE IS AUTOMATICALLY GENERATED BY THIS 
ROUTINE SO THAT *N* POINTS ARE PLOTTED WITH A CONSTANT 
INCREMENT IN X FROM X='LEFT* TO X-'RIGHT'. 

arguments: 

IARRAY * (INTEGER) ARRAY OF Y VALUES TO BE PLOTTED 
N * NO. OF VALUES OF IARRAY TO USE 
(STARTING WITH IARRAY ( 1 ) ) 

IYOFST * (INTEGER) OFFSET (IN TEK. UNITS) TO BE ADDED 
TO EACH Y VALUE AFTER MULTIPLYING BY YSCALE 
YSCALE * (REAL) SCALE FACTOR BY WHICH TO MULTIPLY 
EACH Y VALUE. 

LEFT * (INTEGER) X-COORDINATE FOR LEFT SIDE OF PLOT 
RIGHT * (INTEGER) X-COORDINATE FOR RIGHT SIDE OF PLOT 
BOTH LEFT S RIGHT ARE IN TEKTRONIX SCREEN UNITS 

FOR IYQFST=0 ? YSCALE=1.0? LEFT*0 ? & RIGHT=4095? THE 
VALUES IN IARRAY WILL BE PLOTTED ACROSS THE ENTIRE WIDTH 
OF THE SCREEN? WITH VALUES OF 0 IN IARRAY PLOTTED AT 
THE EXTREME BOTTOM & VALUES OF 3070 PLOTTED AT THE TOP 
OF THE SCREEN. USUALLY IT IS BEST TO AVOID PLOTTING ALL 
THE WAY TO THE EDGES OF THE SCREEN? ESPECIALLY IF HARD 
COPIES ARE DESIRED. 

AUTHOR: WILLIAM G. CROSIER- 

DATE : JULY 1980 

REVISED.* DEC. 1980 

INTEGER I ARRAY < N ) ? IYOFST?LEFT ? RIGHT ? MIN? MAX? IX? IY 
REAL YSCALE ?DELTAX 

DELTAX * FLOAT ( RIGH T-LEFT ) /FLOAT ( N-1 ) !X INCREMENT 

IX = LEFT 

IY = IFIX(FLOAT( IARRAY ( 1 ) )* YSCALE ) + IYOFST 
MOVE TO POSITION OF FIRST POINT 
CALL MPLOT ( I X ? I Y ? 0 ) 

DRAW LINES BETWEEN EACH OF THE N POINTS 
DO 20 1=2 ?N 

IX = LEFT + DELTAX*FLOAT( 1-1 ) 

IY = IFIX(FLOAT ( IARRAY ( I ) ) * YSCALE > + IYOFST 
CALL MPLOT (IX?IY?1) 

CONTINUE 

RETURN 

END 
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to. 





ORIGINAL r.'CZ J 3 

OF POOR QUALITY 




FORTRAN If 

Storase 

Map for Proeram Unit 

ARYPLT 



Local 

Variables? .PSECT 

♦DATA? Size * 000052 

( 21. 

words ) 


Name 

Type 

Offset 

Name Type Offset 

Name 

Type 

Offset 

DEI. TAX 

R*4 

000026 

I 1*2 000032 

IX 

1*2 

000022 

IY 

1*2 

000024 

IYOFST 1*2 0 000004 

LEFT 

1*2 0 

000010 

MAX 

1*2 

000020 

MIN 1*2 000016 

N 

1*2 9 

000002 

RIGHT 

1*2 0 

000012 

YSCALE R*4 <? 000006 





Locsl and COMMON Arrays* 

Name Type Section Offset Si 2 e Dimensions 

I ARRAY 1*2 <? *DATA 000000 **** < #* ) (N) 

Subroutines? Functions? Statement and Pt;ocessor-Def ined Functions* 


Name Type Name 

FLOAT R*4 IFIX 


Type Name Type Name 

1*2 MPL0T 1*2 


Type Name 


Type 
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Subroutine XYPLOT 


This Is a general purpose plotting routine for plotting a real array of Y 
values versus a real array of X values. If the parameter ICOOE Is set to 0, 
then only the points themselves are plotted (with asterisks and If ICOOE 
Is some other value, then the points are connected with straight lines. More 
Information on the arguaents/parameters passed to XYPLOT Is given In the 
program listing. Note that N, L, R, B, T, and ICOOE are all of Integer type, 
but that al 1 other arguments are real . 

An example of the use of XYPLOT, along with GRID and ANOTAT, Is given In 
the third sample program (GRTEST) In another section of this report. Please 
refer there for more Information on a typical use of XYPLOT. 

Routine XYPLOT uses the following subroutine. In addition to those from 
the Fortran Library: 

MPLOT 

The listing for XYPLOT follows. 
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1 


originm. rv -33 
OF poor? Q' ' Y 


FORTRAN IV 


V02. 1-1 


W*d 23-D*c-81 10:23J08 


001 


0001 SUBROUTINE XYPLOT < X t Y » N f L f Rf Bf T» XMIN f MAX f THIN f YMAX f ICCDE > 

ROUTINE TO PLOT AN ARRAY OF Y VALUES VERS'i AN ARRAY OF X VALUES 
ON A TEKTRONIX TERMINAL. THE POINTS ARE CONNECTED WITH STRAIGHT 
LINES. 

arguments: 

X * ARRAY OF X-COORDINAYE VALUES TO PLOT (REAL) 

Y « ARRAY OF Y-COORDINATES (REAL) 

N » NUMBER OF PAIRS OF X-Y VALUES TO PLOT (INTEGER) 

L - LEFT BOUNDARY OF PLOTTING AREA ON SCREEN (INTEGER) 

R » RIGHT PLOT BOUNDARY (INTEGER) 

B * BOTTOM PLOT BOUNDART (INTEGER) 

T - TOP PLOT BOUNDARY (INTEGER) 

XMIN - X-VALUE CORRESPONDING TO LEFT SIDE OF PLOTTING AREA (REAL) 
XMAX » X-VALUE CORRESPONDING TO RIGHT SIDE OF PLOTTING AREA (REAL) 
YMIN - Y-VALUE CORRESPONDING TO BOTTOM SIDE OF PLOTTING AREA (REA 
YMAX » Y-VALUE CORRESPONDING TO TOP SIDE OF PLOTTING AREA (REAL) 
ICODE - CODE FOR CONTROLLING TYPE OF LINES TO DRAW BETWEEN 
POINTS (INTEGER). IF ICODE IS POS THEN LINES WILL BE 
CLIPPED IF THEY WOULD EXTEND PAST THE PLOT BOUNDARIES. 

IF ICODE IS 0 OR NEG THEN THE LINES MAY EXTEND PAST THE 
BOUNDARIES. THE ABSOLUTE VAI.UE OF ICODE DETERMINES THE 
TYPE OF LINES TO DRAW (SEE ROUTINE MPLOT ) • 

IF ICODE IS Of THEN THE POINTS ARE PLOTTED WITH ASTERISKS 
(*)» BUT NO LINES ARE DRAWN BETWEEN THEM. 

examples: 

ICODE=0 WILL PLOT ASTERISKS WITH NO CONECTING LINES. 

IC0DE=1 WILL DRAW NORMAL SOLID LINES BETWEEN POINTS 
2 WILL CLIP LINES TO POINTS OUTSIDE PLOT BOUNDARY. 
IC0DE=-97 WILL DRAW DOTTED LINES BETWEEN POINTS f t 
WILL ALLOW THE LINES TO EXTEND OUTSIDE THE 
DESIGNATED PLOT BOUNDARIES. 

NOTE: ARGUMENTS LfR,Bf I T ARE ALL OF INTEGER TYPE AND ARE 
IN TEKTRONIX SCREEN UNITS. LSR MUST BE BETWEEN 0 AND 4095. 

B2T MUST BE BETWEEN 0 AND 3120. THESE PARAMETERS DETERMINE 
THE PORTION OF THE SCREEN TO BE USED FOR A PLOT. 

ARGUMENTS XMIN f XMAX f YMIN f S YMAX ARE ALL REAL AND ARE IN USER 
UNITS (SAME AS IN ROUTINE ANOTAT). THEY MAY BE IN ANY RANGE 
DESIRED. 

AUTHOR: WILLIAM G. CROSIER 
DATE: FEB. 1981 

c 

'• '>02 INTEGER N f L f R f B f T f ICODE t IX f I Y f IXOFST f I YOFST 

0003 REAL X(N)fY(N)f XMIN f XMAX f YMIN f YMAX f XSCALE f YSCALE 

C 

0004 COMMON /GRFCOM/ MCHRSZ f LS I2E f I WIDTH f IHIGHT f IENHAN f T IMERA f 

l? TIMHDC 


C 

C 

c 

c 

c 

c 

c 

c 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


000 ! 


CALCULATE SCALE FACTORS TO CONVERT X 4 r VAl'JES 
FROM USER UNITS INTO TEKTRONIX SCREEN UNITS 
XSCALE = FLOAT (R-L/ / • XMAX-XMIN) 
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ORIGINAL {3 
OF POOR QUALITY 


FORTRAN TV 


002 , 1-1 


Wed 23-Dec-81 10:23508 


PAGE 002 


0006 

0007 

0008 


IXOFST ■ L - XMIN#XSCALE 
YSCALE = FLOAT (T-EO / (YMAX-YMIN) 
IYOFST * Et - Yh IN* YSCALE 


0009 

0010 
0011 
0012 


IF I CODE IS NOT Of DRAW LINES BETWEEN EACH OF THE N POINTS 
DO 50 1=1 fN 

IX = IFIX (FLOAT < X < I ) ) * XSCALE) + IXOFST 
IY = IFIX<FLQAT<Y<D ) * YSCALE) + IYOFST 
IF (ICODE ,LE. 0) GO TO 30 

PREVENT PLOT FROM EXTENDING PAST DESIRED BOUNDARIES 


0014 


IF <IX .LT. L) IX=L 

0016 


IF (IX .GT, R> IX=R 

0018 


IF < I Y .LT, B) I Y=B 

0020 


IF ( I Y .GT. T) I Y=T 

0022 

30 

CONTINUE 


C IF 

ICODE=0. PLOT POINTS ONLY WITHOUT CONNECTING 

0023 


IF (ICODE .EQ. 0) GO TO 40 

0025 


I TYPE = 0 

0026 


IF (I.GT.l) ITYPE = IABS< ICODE) 

0028 


CALL MPLOTdX, IYf ITYPE) 

0029 


GO TO 50 

0030 

40 

CALL HPLOT < IX-I WIDTH/2 ? IY-IHIGHT/3f-l ) 

0031 


TYPE 45 

0032 

45 

FORMAT ('+*') 

0033 

50 

CONTINUE 

0034 


RETURN 

0035 


END 


LINES 
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ORKSliWU. ? .. .f;7 
OF POOR QUALITY 


FORTRAN IV 


Stors.e 

Hap for ProSram Unit XYPLOT 



Local 

Variables* .PSECT 

♦DATA* Size = 

> 000100 < 32. words) 


Name 

T ype 

Offset 

Name Type 

Offset Name 

T ype 

Offset 

B 

1*2 

II 

000012 

I 1*2 

000054 I CODE 

1*2 <? 

000026 

I TYPE 

1*2 


000056 

IX 1*2 

000034 IXOFST 

1*2 

000040 

IY 

1*2 


000036 

IYOFST 1*2 

000042 L 

1*2 0 

000006 

N 

1*2 

<? 

000004 

R 1* 2 <? 

000010 T 

1*2 Q 

000014 

XMAX 

R*4 

e 

000020 

XhIN R*4 0 

000016 XSCALE 

R*4 

000044 

YHAX 

R*4 

e 

000024 

YMIN R*4 (1 

000022 YSCALE 

R*4 

000050 


COHHON 

Block 

/Gf'FCOH/ * 

Size * 000022 

< 9 * 

words ) 



Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

MCHRSZ 

1*2 

000000 

LSIZE 

1*2 

000002 

I WIDTH 

1*2 

000004 

IHIGHT 

TIHHDC 

Local , 

1*2 000006 
R*4 000016 

and COMHON Arrays 

IENHAN 

• 

♦ 

1*2 

000010 

T I HERA 

R*4 

000012 

Name 

X 

Y 

Type Section 

R*4 @ «DATA 

R*4 (1 SDATA 

Offset 

000000 

000002 

**** 

**** 

-Size 

( ** 

( ** 

— Dimens 
) < N ) 

) < N ) 

ions 



Subroutines* Functions* Statement and Processor-Def ined Functions* 

Name Type Name Type Name Type Name Type Name Type 

FLOAT R*4 I ABS 1*2 IFIX 1*2 HPLOT 1*2 
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Subroutine GRID 


This subroutine Is used to draw horizontal and/or vertical grid lines 
over a specified portion of the terminal screen. All of the arguments/ 
parameters are of Integer type, and their functions should be clear from 
reading the program listing below and the third example program (GRTEST) In 
another section of this report. The values for arguments LEFT, RIGHT, BOTTOM , 
and TOP should normally be the same In the calls to XYPLOT (or ARYPLT), GRID, 
and ANOTAT. If NXDIV Is set equal to 0, then no vertical grid lines will be 
drawn. Similarly, If MYDIV Is set equal to 0, then no horizontal grid lines 
will be drawn. ICODE specifies the line type, as described In the section on 
routine MPLOT. 

Routine GRID uses the following subroutine. In addition to those from the 
Fortran Library: 

MPLOT 

The program listing for GRID follows. 
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FORTRAN IV 


V02. 1-1 


PAGE 001 


ORIGIN . 

OF P -'" r 


Wed 23-Dec~81 10 .*23: 21 


0001 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

0002 

0003 

c 

c 

c 

0004 

0006 

0007 

0008 

0009 

0010 

0011 

0012 50 
c 
c 


SUBROUTINE GRID( NXDIV, NYDIV * LEFT , RIGHT * BOTTOM * TOP * ICODE > 
THIS ROUTINE DRAWS GRID LINES OVER A DESIRED PORTION 
OF THE SCREEN ON A TEKTRONIX TERMINAI . 

ARGUMENTS (ALL ARE INTEGERS): 

NXDIUss# OF X-AXIS DIVISIONS = # OF VERT. LINES - 1 
NYDIV=* OF Y-AXIS DIVISIONS = * OF H0RI2. LINES - 1 
LEFT = LEFT BOUNDARY OF PLOTTING AREA IN TEK. 

SCREEN UNITS (MIN, X VALUE) 

RIGHT = RIGHT BOUNDARY (MAX. X) 

BOTTOM = LOWER BOUNDARY (MIN. Y) 

TOP = UPPER BOUNDARY (MAX Y) 

ICODE = DESIGNATES TYPE OF LINES TO DRAW 

(SEE ROUTINE MPLOT OR TEKTRONIX MANUAL) 

EX. : 

ICODE = 1 — DRAW SOLID LINES 
ICODE = 97 — DRAW DOTTED LINES 

author: william g. crosier 

date: 9 DEC. 1980 

INTEGER NXDIV* NYDIV * LEFT* RIGHT * BOTTOM* TOP* ICODE* IX * I Y 
REAL DELTA 

DRAW HORIZ. LINES 

IF (NYDIV .LE. 0) GO TO 80 
CALL MPLOf (LEFT *B0TT0M*O) 

CALL MPLOT ( RIGHT * BOTTOM* ICODE ) 

DELTA=FLOAT ( TOP-BOTTOM ) /FLOAT ( NYDIV ) 

DO 50 K® 1 * NYDIV 

IY = BOTTOM + IFIX (DELTA * FLOAT (K)) 

CALL MPLOT (LEFT * IY*0) 

CALL MPLOT (RIGHT* IY* ICODE) 

DRAW VERTICAL LINES 


0013 80 

0015 

0016 

0017 

0018 

0019 

0020 

0021 100 
0022 999 

002 3 

C 


IF (NXDIV .LE. 0) GO TO 999 
CALL MPLOT (LEFT* BOTTOM *0) 

CALL MPLOT (LEFT, TOP* ICODE) 

DELTA = FLOAT (RIGHT-LEFT) / FLOAT (NXDIV) 
DJ 100 K=1*NXDIV 

[X = LEFT + IFIX(DELTA*FLOAT(K) ) 

CALL MPLOT ( IX* B0TT0M*0) 

CALL MPLOT (IX* TOP* ICODE) 

RETURN 

END 
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ORIGINAL 

OF POOR QUALSHY 


FORTRAN IV 

Storage 

Map for 

Program Unit GRID 



l.ocsl 

Variables* .PSECT 

♦DATA* 

Size ® 000044 ( 

18* words ) 


Name 

Type 

Offset 

Name 

Type Offset 

Name 

Type 

Offset 

BOTTOM 

1*2 e 

000010 

DELTA 

R*4 000022 

I CODE 

1*2 (? 

000014 

IX 

1*2 

000016 

IY 

1*2 000020 

K 

1*2 

000026 

LEFT 

1*2 Q 

000004 

NXDIV 

1*2 1? 000000 

NYDIV 

1*2 l? 

000002 

RIGHT 

1*2 <? 

000006 

TOP 

1*2 <? 000012 




Suhrout ines f 

Functions * 

Statement and Processor- 

-Defined 

1 Functions? 

Name 

Type 

Name Type Name Type Name 

T ype 

Name 

Type 

FLOAT 

R*4 

IFIX 1*2 MPLOT 1*2 
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Subroutine ANOTAT 


This subroutine Is used to anotate (label) the horizontal and/or vertical 
plot axes with nunerlcal user units at some or all grid lines (previously 
produced by routine GRID). Most of the arguments or parameters for this 
routine are the same as for GRID and XYPLOT, and you may refer to the listings 
of those routines and of ANOTAT Itself (given below) for more Information. In 
addition, the third sample main program (GRTEST) In another section of this 
report gives an example of how these routines can be used. You may note that 
the values for NXDIV and NYDIV may be different In the call to ANOTAT from 
what they were In the call to GRID. This can be done If you want to label 
only every second, fifth, etc., grid line with ANOTAT. 

Be sure to call routine CHRSIZ at least once In your program before 
calling ANOTAT. This Is necessary so that ANOTAT will know what the current 
character size Is, so that It can position the nunerlcal units for the axes 
properly. In addition, make sure that arguments L and B are large enough so 
that ANOTAT will not attempt to type the nunerlcal units to the left of, or 
below, the allowable plotting area. This may especially be a problem If one 
of the two larger character sizes are being used. If L or B are too small, 
then the numbers may be typed over the plot axes and be difficult to read. 

ANOTAT requires the following subroutine, In addition to those from the 
Fortran Library: 

MPLOT 

A listing of ANOTAT follows. 
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FORTRAN IV 


V02. 1-1 


PAGE 001 


0001 


0002 

0003 

0004 


0005 


0006 

0003 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

00 1 * 


Oh;- 

OF POO l O '- 


Wed 23-Dec-81 10! 23! 32 

SUBROUTINE ANQTAT ( NXDIV f NYDIV f L f R f B f T f XMIN f XHAXf YMIN f YMAX ) 

C 

C ROUTINE TO ANOTATE PLOT AXES WITH USER UNITS AT GRID LINES. 

C IF ROUTINE GRID IS ALSO USED t THE FIRST 6 ARGUMENTS SHOULD BE 
C IDENTICAL TO THOSE USED IN GRID. IN THAT CASE* THE GRID LABELS 
C WILL BE TYPED JUST OUTSIDE THE PLOTTING AREA. 

C 

C ARGUMENTS! 

C FIRST 6— SAME AS FOR GRID (ALL INTEGER)! 

C NXDIV=# X-AXIS DIVISIONS = * VERT. LINES - 1 

C NYDIV=# Y-AXIS DIVISIONS = # HORIZ LINES - 1 

C L = LEFT BOUNDARY OF PLOT IN TEKTRONIX SCREEN UNITS 

C R = RIGHT BOUNDARY 

C B = BOTTOM BOUNDARY 

C T = TOP BOUNDARY 


c 

XMIN = 

MIN. 

X 

VALUE 

TO 

BE 

TYPED 

AT 

BOTTOM 

LEFT 

(REAL) 

c 

XMAX = 

MAX. 

X 

VALUE 

TO 

BE 

TYPED 

AT 

BOTTOM 

RIGHT 

(REAL) 

c 

YMIN = 

MIN. 

Y 

VALUE 

TO 

BE 

TYPED 

AT 

BOTTOM 

LEFT 

(REAL) 

c 

YMAX = 

MAX. 

Y 

VALUE 

TO 

BE 

TYPED 

AT 

TOP LEFT (REAL) 


C 

C IF NXDIV IS GREATER THAN Of THEN NXDIV+1 NUMBERS (STARTING WITH 
C XMIN S ENDING WITH XMAX) ARE TYPED BELOW VERT. GRID LINES (DRAWN 
C SEPARATELY WITH ROUTINE GRID) WITH AN F6.D FORMAT JUST BELOW PLOT. 
C IF NYDIV IS GREATER THAN Of THEN NYDIV+1 NUMBERS (STARTING WITH 
C YMIN g ENDING WITH YMAX ) ARE TYPED TO THE LEFT OF HORIZ. GRID 
C LINES (DRAWN WITH GRID) WITH AN F6.D FORMAT. 

C 

C AN F6.0 FORMAT WILL BE USED UNLESS THE MIN. X OR Y VALUE IS 
C GREATER THAN -9.99 AND THE MAX VALUE IS LESS THAN 99. 9f IN 
C WHICH CASE AN F6.3 FORMAT WILL BE USED INSTEAD. 

C 

C AUTHOR! WILLIAM G. CROSIER 
C DATE! 16 DEC. 1980 

C 

INTEGER NXDIVfNYDIVfLfRfBfTf IX f IY 

REAL XMIN f XMAX f YMIN f YMAX t DELTA , DELTA1 fXfYf FMT ( 3 ) 

COMMON /GRFCOM/ MCHRSZ f LSIZE f I WIDTH f IHIGHT f IENHAN f 

e timerafTimhdc 

c 

DATA FMT /4H('+' f4HfF6. f4H0> / 

C 

C LABEL Y (VERTICAL) AXIS 

C 

IF (NYDIV ,LE, 0) GO TO 80 
FMT ( 3 ) = -0) ' 

IX = L - 6*1 WIDTH 

DELTA = FLOAT ( T-B ) / FLOAT (NYDIV) 

DELTA1 = (YMAX-YMIN) / FLOAT (NYDIV ) 

DO 50 K=1fNYDIV+1 

IY = B + IFIX( DELTAKFLQAT ( N-l ) > - IFIX(0. 2# IHIGHT) 

C MOVE TO DESIRED POSITION TO LEFT OF AXIS. ALPHA MODE 

CALL MPLOTIIXf IYf-1) 

Y = YMIN + DEL IA1 * FLOAT ( K-l > 

IF 1 YMIN , GT .-9,99 .AND. Y MAX. GT. -9. 99 .AND. rMIN . L T . 99 . 9 
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•AND* YMAX.LT.99.9) FMT < 3 >* ' 3 ) 


page oo: 


0018 


TYPE FMT, Y 

0019 

50 

C 

c 

c 

80 

CONTINUE 


LABEL X (HORIZ.) AXIS 

0020 

IF (NXDIV .LE. 0) GO TO 999 

0022 


FMT ( 3 ) * '0) ' 

0023 


IY = B - 1 • 2*IHIGHT 

0024 


DELTA = FLOAT < R-L ) / FLOAT (NXDIV ) 

0025 


DELTA1 » (XMAX-XMIN) / FLOAT (NXDIV) 

0026 


DO 100 K=1,NXDIV+1 

0027 


IX ■ L + IFIX < DELTA*FLQAT (K-l ) ) - 4*IWIDTH 


C 

MOVE TO DESIRED POSITION BELJ3W AXIS, ALPHA i 

0028 


CALL MPLOT <IX»IY»-1) 

0029 


X = XMIN + DELTA1 * FLOAT < K-l ) 

0030 


IF ( XMIN ♦ GT . -9 . 99 .AND* XMAX . GT . -9 . 99 .AND, 
0 .AND. XMAX. LT. 99. 9) FMT<3)='3 

0032 


TYPE FMT, X 

0033 

100 

P 

CONTINUE 

0034 

U 

999 

RETURN 

0035 


END 


C 


FORTRAN IV Storage Map for Program Unit ANOTAT 


Local 

Variables, .PSECT 

♦DATA, 

Size ■ 

= 000132 ( 

45. words) 


Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

B 

1*2 

0 

000010 

DELTA 

R*4 

000044 

DELTA1 

R*4 

000050 

IX 

1*2 


000040 

IY 

1*2 

000042 

K 

1*2 

000064 

l 

1*2 

0 

000004 

NXDIV 

1*2 0 

000000 

NYDIV 

1*2 0 

000002 

R 

1*2 

0 

000006 

T 

1*2 0 

000012 

X 

R*4 

000054 

XMAX 

R*4 

0 

000016 

XMIN 

R*4 0 

000014 

Y 

R*4 

000060 

YMAX 

R*4 

0 

000022 

YMIN 

R*4 0 

000020 





COMMON 

Block 

/GRFCOM/ , 

Size = 

000022 

( 9. 

words ) 



Name 

Type 

Off set 

Name 

Type 

Offset 

Name 

T ape 

Offset 

MCHRSZ 

1*2 

000000 

LSIZE 

1*2 

000002 

IWIDTH 

1*2 

000004 

IHIGHT 

TIMHDC 

1*2 

R*4 

000006 

000016 

IENHAN 

1*2 

000010 

TIMERA 

R*4 

000012 


Local and COMMON Arrays l 

Name Type Section Offset Size Dimensions 

^MT R*4 $DATA 000024 000014 < 6.) (3) 


Subroutines* Functions* Statement and Processoi — Defined Functions J 

Name Type Name Type Name Tyre Name Type Name Type 

FLOAT f\'*4 IF IX 1*2 MPLOT 1*2 
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Subroutine SELL 


This subroutine Is used to make the terminal beep or ring Its bell. The 
sound produced depends on the particular terminal. The duration of the sound, 
as well as Its modulation, can be controlled with the two parameters or 
arguments NUMBER and IDELAY. Examples of the use of this routine are given In 
the third sample main program ( GRTEST ) , In another section of this report. 

Routine BELL requires the following subroutines, In addition to those 
from the Fortran Library: 

WAIT 

ITTOUR (from the System Subroutine Library) 

The listing for subroutine BELL follows. 
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0001 


0002 

0003 

0005 

0006 
0007 

0009 

0010 
0011 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


40 

100 


c 


SUBROUTINE BELL ( NUMBER? I DELAY ) 

RING TERMINAL BELL/BEEP WITH VARIABLE DURATION 
1 MODULATION CONTROL 

NUMBER » NO. OF BELL CHARACTERS TO TRANSMIT 
(CONTROLS DURATION) 

IDELAY « NO. OF 1/60 SEC INCREMENTS TO 

WAIT BETWEEN BELLS (CONTROLS MODULATION 
1 PERCEIVED FREQUENCY). 

IDELAY CAN BE 0 FOR A CONTINUOUS TONE. WITH DURATION 
CONTROLLED BY NUMBER. OR IDELAY CAN BE A POSITIVE 
INTEGER TO PRODUCE A BUZZING SOUND OR DISCRETE BEEPS. 

NOTE: SOUND IS DEPENDENT ON THE TERMINAL * ON ITS 
BAUD RATE SETTING. 

DO 100 K =1. NUMBER 

IF (IDELAY . LT . 1) GO TO 40 
DELAY * FLOAT( IDELAY) / 60.0 
CALL WAIT(DELAYfO) 

IF ( IT TOUR ( 7 ) .NE. 0) GO TO 40 ! SEND BELL 

CONTINUE 
RETURN 
END 


FORTRAN IV Storaae Map for Proaram Unit BELL 


Local 

Variables. .PSECT 

♦DATA. 

Size = 000012 ( 

5 • 

words ) 


Name 

Type Offset 

Name 

Type Offset 

Name 

Type 

Off set 

DELAY 

NUMBER 

R*4 000006 

1*2 <? 000000 

IDELAY 

1*2 <? 000002 

N 

1*2 

000004 


Subroutines. Functions. Statement and Processor-Defined Functions! 


Name 

T ype 

Name 

T ype 

Name 

T ype 

Name 

T ype 

Name 

FLOAT 

R*4 

IT TOUR 

1*2 

WAIT 

R*4 
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Subroutine WAIT 


This routine uses the RT-ll system line frequency clock to time a waiting 
period. The u»er simply passes the routine a real value In argument/ parameter 

SEC which specifies the duration, In seconds, of the waiting period. The 

normal functions of the line time clock are not affected. When the waiting 
period has elapsed, then control Is retu. .ied to the calling program. If you 
want to be able to terminate the wait prematurely (In less than "SEC" 
seconds), then an Interrupt service routine can be used to set the argunent 
I ABORT equal to a non- zero value. You may want to do this In a real-time 
experiment control program. If something happened during a programmed wait or 
if you pressed a button on a control panel, for example. If you do not need 
to prematurely terminate the waiting period, then set IAB0RT equal to 0. An 
example of the use of this routine can be seen In the listing for subroutine 
ERASE. This routine calls WAIT after sending the command to erase the 

terminal screen, so that the terminal will have enough time (normally TIM ERA 
Is 1.5 seconds) to completely clear the screen. 

Subroutine WAIT may also be used to time the periods between data 

acquisition samples. If the sampling rate If fairly slow (60 HZ or slower). 
The accuracy and resolution of the programmed wait Is 1/60 > 0.017 second, 
since that Is the time between cycles of the line frequency. There should be 
no cumulative time error between waiting periods, however, as long as not too 
many CPU operations are performed, so that the long term accuracy of the times 
measured should be quite good. As an example of using WAIT to control data 
acquisition, the following could be used to get 1000 samples of signals on 
analog channel 3 and 5 with the samples taken every 1/30 ■ 0.033 second: 

INTEGER IDATA(IOOO) ,IDATB(1000) 

IPGNCD * 2 ! PROGRAMMABLE GAIN COOE 

00 10 K- 1,1000 IDO 1000 TIMES 

IDATA(K) « ISAMPA(3 ,IPGNCD,0) ! SAMPLE CHANNEL 3 

IDATB(K) * ISAMPA( 5, IPGNCD, 0) I SAMPLE CHANNEL 5 

CALL WAIT (0.033,0) 

10 CONTINUE 
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At the end of 33 seconds, 1000 samples would be collected from Channel 3 
and stored In array I DATA, while 1000 samples from Channel 5 would be In array 
IDATB. If numerous computations such as averaging of many samples are per- 
formed between waiting periods, then the period may actually be longer than 
desired, since the computations make take more than 1/60 of a second. 

If you have a 50 HZ line frequency system clock, rather than the 60 HZ 
usually used In the United States, be sure to change this line In the program: 

from: TICKS » SEC*60. + 0.5 
to: TICKS ■ SEC*50 . ♦ 0.5 

Oo this only If you have a 50 HZ line frequency clock. 

Routine WAIT requires the following subroutines from the System 
Subroutine Library: 

GTIM (Time of Day In clock ticks past midnight) 

JADD (IntegarM addition) 

JAFIX (Real *4 to Integer*4 conversion) 

JCMP (IntegerM compare) 

JJCVT (Interchange halves of Integer*4 variable) 


A listing of WAIT follows 


ORIGINAL PAG* IS 
OF POOR QUALITY 

FORTRAN IV V02.1-1 Wed 23-Dec-8l 10: 23 MS PAGE 001 


0001 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

0002 

0003 

0004 

0005 

0006 

0007 

0008 10 

0009 

0010 
0012 

0014 99 

0015 

C 


SUBROUTINE WAIT < SEC » I ABORT ) 

LINE TIME CLOCK WAIT ROUTINE 
WRITTEN BY: WILLIAM G. CROSIER 

SEC * NUMBER OF SECONDS (REAL. NOT INTEGER) TO WAIT 
USES LINE TIME CLOCK FOR TIMING CONTROL. 

THE J SUBROUTINES USED HERE PERFORM INTEGER*4 ARITHMETIC. 

RESOLUTION * ACCURACY * APPROX. 0.017 * i/*0 SECOND 
THIS ROUTINE DOES NOT AFFECT NORMAL FUNCTIONS OF LTC . 

IF PARAMETER I ABORT BECOMES NON-ZERO DURING 

THE WAITING PERIOD ( IF SET BY AN INTERRUPT ROUTINE) » 

THEN THE WAIT IS IMMEDIATELY TERMINATED. 

IF THIS FEATURE IS NOT NEEDED » USE A VALUE OF 0 FOR IABORT 
INTEGERS ITIM1.ITIM2.IDELTA 

CALL GT [M < I TIM1 > 'STORE CURRENT TIME IN I TIME 

CALL JJCVT(ITIMl) ! INTERCHANGE WORDS 

TICKS * SEC#60. + 0.5 'CONVERT SEC TO CLOCK TICKS 

CALL JAFIX ( TICKS » I DELTA ) ! CONVERT TO INTEGER- 

CALL JADD< ITIril » I DELTA* ITIM1 ) ‘CALCULATE STOP TIME 

CALL 6TIM( ITIM2) > GET CURRENT TIME OF DAY 

CALL JJCVT < ITIM2 ) 'INTERCHANGE WORDS 

IF (IABORT .NE. 0) GO TO 99 (CHECK FOR ABORT 

IF ( JCMP(ITIM2f ITIM1 ) .LT. 0) GO TO 10 

RETURN 

END 


FORTRAN IV S to 

Local Variables. .P 

Name Type Offset 
IABORT I *2 8 000002 
ITIM2 1*4 000010 


e Map tor Program Unit 

T 4DATA. Size - 000024 

Name Type Offset 
IDELTA 1*4 000014 

SEC R*4 l? 000000 


WAIT 

( 10. words) 

Name Type Offset 
ITIM1 1*4 000004 

TICKS R*4 000020 


Subroutines * 


Functions. Statement and Froces sor-Def ineo Functions: 


Name Type 
Oflrt R*4 


Name Type 
JADD 1*2 


Name Type 
JAFIX 1*2 


Name Type 
JCMF' 1*2 


Name Tyre 
JJCVT 1*2 
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Subroutine I SAMP A 


This Is a Fortran-callable subroutine, written In MACRO Assembly 
language, for sampling an analog signal with an analog-to-dlgltal (A/D) 
converter. It was written In MACRO so that It can execute as quickly as 
possible, but still be usable with FORTRAN programs. With this routine, 
sampling rates of several hundred samples per second can easily be achieved, 
even with some computations performed between samples. For accurate control 
of the time Intervals between samples, you may use routine WAIT If the 
sampling rate Is 60 HZ or slower. Otherwise, you should use a programmable 
clock/timer such as the KW-11P. If ISAMPA Is called by a MACRO Interrupt 
service routine for the KW-11P, make sure that the routine uses the normal 
PDP-11 Fortran calling conventions for passing arguments, etc. 

If accurate control of tho sampling rate Is not a requirement, but you 
need to sample a large number of values In a certain time period and average 
them In order to reduce noise effects, you may use the following procedure. 
First, get the current time of day (In seconds past midnight) with the RT-11 
system routine SECNDS, or wait until an appropriate external event occurs. 
Second, call ISAMPA, convert the returned sampled value to real or double 
precision, and add it to a real or double precision variable used as an 
accumulator. Repeat the sampling and accumulating untl'.’ either you have 
enough samples, or until enough time has elapsed. (Use the SECNDS subroutine 
again.) Finally, divide by the number of samples collected. A real variable 
(rather than integer) should be used If you are using a 12-bit A/D converter 
and are adding together more than 16 samples, because a 16-bit integer 
accumulator can be overflowed by adding more than 16 12-bit values together If 
each of them are near full scale, generally, no error message will occur if 
this happens, since PDP-11 Fortran does not check for an overflow on an 
integer add operation. S1m1 ; ar1y, you should use a double precision 
accunulator If you add together more than about 2000 samples, because you can 
drop bits when doing so with Real *4 arithmetic. 

For examples of how ISAMPA may be used, refer to the discussions for 
subroutine WAIT and for the first sample main program (ADTEST), In another 
section of this report. 
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If you are using a DEC ADV-11A A/D converter, then you must mask out the 
four most significant bits, since DEC uses them for other purposes. In 
addition, DEC'S ADV-11A converters can only be used with an offset binary 
format, so that a value of 4000 (octal) or 2048 (decimal) must be subtracted 
from the sampled value In order to convert It to the normal two's complement 
coding. The following will mask out the 4 MSB's and convert the value to 2's 
compl ement: 

I-C7777 .AND. ISAMPA(ICHAN,0,D) •.“4000 

This Is necessary only with DEC A/D boards. Note also. In the above 
example, that a value of 0 must be used for the second argument ( IPGNCD) , 
since the DEC boards do not have programmable gain. 

ISAMPA requires no subroutines. 

A listing of ISAMPA follows. 
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hIiCSK 1-34# 1-35 

AODrlTA 1-35# 1-51 

Utl'EC 1-42 1-47# 

UTIiAT 1-46 1-51* 

JSAnPA 1-33 1-37# 


Subroutine DISK 10 


This routine Is used to read or write binary unformatted data. In a 
random fashion, to a sequential disk file. Variable length records are 
supported (In multiple of 256 Integer words), and the routine Is somewhat more 
economical of memory and CPU time than the normal Fortran disk I/O, since the 
transfer takes place directly from the arrays In the user's program, rather 
than through Intermediate buffers. Any binary data can be transferred, 
regardless of whether the calling program treats It as logical. Integer, real, 
or string data. The principal restrictions on the data Is that It must all be 
In contiguous memory locations, and that only multiples of 256 words should be 
transferred normally through each call to DISKIO. Generally, the data should 
all be placed In a COMMON block In order to force the compiler to place It all 
In contiguous locations, unless It Is In a single array. 

The arguments/parameters for this routine are discussed In the program 
Using below. The file name passed In argument FILNAM can be any valid RT-11 
file name. Note, however, that It must be exactly 12 characters (bytes) long, 
with no colons or periods within It to separate the device name or file type/ 
suffix. Trailing spaces are allowed, however, at the end of each portion of 
the file name in order to make the device name identifier exactly 3 characters 
long and the main part of the file name exactly 6 characters long. A null or 
zero byte should follow the 12 character file name (as In the MACRO assembler 
ASCIZ construction). 

Data written by this routine can only be read back with the same routine, 
and not by Fortran READ statements. In addition, this subroutine was designed 
to work only with the RT-11 operating system. 

For more information and examples of the use of DISKIO, refer to the 
fourth sample main program, DISKRW, in another section of this dociment. 

Refer to the program listing which follows for more Information. 
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ORIGINAL r/'"' 13 
OF POOR QUALITY 

FORTRAN IV V02.1-1 Wed 23-Dec-8l 10:23.‘55 PAGE 001 

0001 SUBROUTINE DISKIO ( FILNAM . MODE . BUFFER , NWRDS . IBLK.NBLK. I ERR) 

purpose: READ OR write binary data to a disk file 

WRITTEN BY: WILLIAN 6. CROSIER 
DATE: 11 JUNE 1980 

arguments: 

FILNAn£=ARRAY CONTAINING ASCII FILE NAME (12 CHAR.) 

(FILNAM IS IGNORED WHEN MODE WAS NEG. ON THE LAST 
CALL TO DISKIO.) 

EX: DX1TEST07DAT 

DK MYFILE 
DK FILE7 

M0DE=+1 OR -1 TO CREATE A NEW FILE 8 WRITE OUT DATA TO IT 

(+1 OR -1 WILL CAUSE ANY FILE WITH THE SAME NAME WHICH 
PREVIOUSLY EXISTED TO BE DELETED WHEN THE NEW 
FILE IS CLOSED) 

=2 OR -2 TO MODIFY AN EXISTING FILE (OVERWRITE ALL OR PART) 
=3 OR -3 TO READ DATA FROM AN EXISTING FILE 
IF MODE IS POSITIVE. THE FILE IS CLOSED AFTER THE I/O. 

IF MODE IS NEGATIVE. THE FILE IS NOT CLOSED. SO THAT THE NEXT 
CALL TO DISKIO WILL NOT REQUIRE RE-OPENING THE FILE. 

(THE NEXT CALL WILL ALSO IGNORE FILNAM SINCE THE 
PREVIOUSLY SPECIFIED NAME WILL BE USED AGAIN.) 

IF AN ERROR OCCURS WHEN MODE IS NEG,. THE NEXT I/O 
MAY NOT BE VALID SINCE THE FILE MAY NOT BE OPENED 
PROPERLY, TO AVOID THIS PROBLEM. DO A READ OPERATION 
WITH M0DE=3 TO CLOSE THE FILE IF AN ERROR OCCURS 
WHEN MODE IS NEG. 

note: all FILES ARE unconditionally closed when 

THE PROGRAM TERMINATES. REGARDLESS OF WHETHER MODE 
WAS POS. OR NEG. ON THE LAST CALL. 

BUFFER=AREA IN MEMORY WHERE DATA IS TO BC TRANSFERRED TO/FROM 

NWRDS=NO. OF INTEGER WORDS TO READ FROM OR WRITTEN INTO BUFFER 
(SHOULD BE A MULTIPLE OF 256) 

IBLK=STARTING BLOCK NO. IN FILE WHERE DATA TRANSEER IS TO OCCUR 

NBLK=NO. OF 256-WORD INTEGER BLOCKS TO ALLOCATE FOR A NEW FILE 
( NBLK IS IGNORED EXCEPT WHEN M0DE=1 OR -1) 

IERR=ERROR CODE RETURNED BY DISKIO 
=0 MEANS NO ERRORS OCCURRED 
= 1 MEANS QUEUE E! EmENT FAILURE OCCURRFD 
=2 MEANS NO I/O CHANNEL WAS AVAILABLE 

-3 MEANS HANDLER FOR SPECIFIED DEVICE CAN'T BE LOADED 
=4 nEANS FILE ALLOCATION FAILED WHEN CREATING FILE 
=5 MEANS A DATA OUTPUT ERROR OCCURRED 
=6 MEANS A FILE LOOKUP FAILURE (COULD NOT FIND FILE; 
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ORIGINAL P*-.£ fcJ 
OF POO*> O' ‘A’ iTY 

FORTRAN 10 002.1-1 Wed 23-Dec-8t 10:23155 PAGE 002 


0002 

0003 

0004 

0005 

0006 
0007 

0009 

0010 


0012 

0014 

0015 

0016 
0017 

0019 

0020 
0021 

0023 

0024 

0025 

0027 

0029 

0030 

0031 

0033 

0035 

0036 


0037 

0039 

0040 

0041 

0043 

0044 

0046 

0047 

0048 

0049 

0050 


C =7 MEANS A DATA INPUT ERROR OCCURRED 

C 

C THIS ROUTINE TAKES CARE OF OPENING * CLOSING THE FILE FOR 

C EACH DATA TRANSFER . SETTING QUEUE ELEMENTS APPROPRIATELY. 

C GETTING AN I/O CHANNEL . FETCHING THE DEVICE HANDLER . CREATING 

C THE FILE ENTRY. 3 DOING THE ACTUAL DATA I/O 

C 

BYTE FILNAM< 12 ) 

INTEGER MODE » BUFFER » NWRDS . I BLK » NBLK . IFI LE < 4 > » FLAG * I ERR » I CHAN » 

8 IPMODE 

COMMON /DISCOM/ IPMODE 
DATA IPMODE /O/ 

IERR=0 

IF < IPMODE.LT. 0) GO TO 60 ! FILE LEFT OPEN? 

C CONVERT FILE NAME TO RADIX-50 

CALL IRAD50 ( 12 > FILNAM. IFILE ) 

IF (IPMODE .NE. 0) GO TO 20 

C FIRST TIME ROUTINE HAS BEEN CALLED » SO SET QUEUE ELEMENTS 
IF ( IQSET ( 2 ) .EQ. 0) GO TO 20 

IERR=1 ! ERROR-QUEUE ELEMENT FAILURE 

GO TO 999 

20 ICHAN=IGETC< IDUMMY ) ! GET AN I/O CHANNEL 

IF (ICHAN . GE . 0) GO TO 30 

IERR=2 ! ERR-NO CHAN. AVAIL, 

GO TO 99 

30 IF <IFETCH<IFILE<1 ) ) , EQ. 0) GO TO 40 ! FETCH DEVICE HANDLE; 

IERR=3 ! ERR-CANNOT LOAD HANDLER 

GO TO 90 

40 IF ( IABS < MODE ) . GT . 1) 60 TO 50 
C CREATE NEW FILE ENTRY 

IF (IENTERUCHAN.IFILE.NBLK) ,GE.O> GO TO 60 

IERR=4 ! ERR-FILE ALLOCATION FAILED 

GO TO 90 

60 IF ( IABS < MODE ) .EQ. 3) GO TO 70 
C WRITE OUT DATA FROM BUFFER 

IF < I WRI TW ( NWRDS » BUFFER > IBLK » ICHAN ) . GE . 0 ) GO TO 90 
IERR=5 ! ERR-DATA OUTPUT 

GO TO 90 

C FIND EXISTING FILE 

50 IF <LOOKUP< ICHAN. IFILE) .GE. 0) GO TO 60 

IERR=6 ! ERR IN FILE LOOKUP 

GO TO 90 

C READ DATA INTO BUFFER- 

70 IF< IREADWCNURDS. BUFFER. IBLK. ICHAN) .GE.O) GO TO 90 

IERR-7 ! ERR IN READING DATA 

90 IF (MODE .LT. 0) GO TO 99 ! LEAVE CHAN. OPEN? 

CALL CLOSEC< ICHAN) 'CLOSE THE I/O CHANNEL 

CALL IFREEC< ICHAN) & FREE IT 

99 IPMODE*MQDE 

999 RETURN 
END 
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OF POOR QUALu / 


FORTRAN 

1 IV 

Storage 

Map for Program Unit 

DISK 10 



Local Variables* .PSECT 

5DATA * 

Size ■ 

( 000034 

( 14. 

woros ■ 


Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Of f set 

BUFFER 

1*2 0 

000004 

FLAG 

1*2 

000026 

IBLK 

1*2 0 

000010 

ICHAN 

1*2 

000030 

I DUMMY 

1*2 

000032 

IEPP 

1*2 0 

000014 

MODE 

1*2 0 

000002 

NBLK 

1*2 0 

000012 

NURDS 

1*2 0 

000006 

COMMON 

Block 

/DISCOM/. 

Size * 

000002 

( 1 . 

words) 



Name 

Type 

Offset 

Name 

Type 

Offset 

Name 

Type 

Offset 

IPMODE 

1*2 

000000 








Local and COMMON Arrays! . 


Name Type Section Offset Size Dimensions 

FlLNAM L*1 0 *DATA 000000 000014 ( 6.) (12) 

IFILE 1*2 SDATA 000016 000010 < 4.) (4) 


Subroutines t Functions* Statement and Processor-Defined Functions! 


Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

Name 

Type 

CLOSEC 

R*4 

IABS 

1*2 

IENTER 

1*2 

IFETCH 

1*2 

IFREEC 

1*2 

ICETC 

1*2 

IQSET 

1*2 

IRAD50 

1*2 

IREADU 

1*2 

IWRITW 

1*2 


LOOKUP 1*2 
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